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PREFACE 



This publication is intended to provide a reference 
to the language of 7090/7094 FORTRAN IV (IB FTC), 
which is a part of the 7090/7094 IBJOB Processor. 
The other parts of the IBJOB Processor are de- 
scribed in IBM 7090/7094 Programming Systems: 
IBJOB Processor , Form C28-6275, and MAP (Mac- 
ro Assembly Program) Language , Form C28-6311. "^. ) 

7090/7094 FORTRAN IV (hereafter referred to as \ 4^ J 

FORTRAN) is a language component of IBJOB. The 
basic concepts of a FORTRAN language are dis- 
cussed in the FORTRAN general information manual, 
Form F28-8074-1. The level of information on the 
following pages is such that the reader must be 
thoroughly familiar with the information in the FOR- 
TRAN general information manual. 

FORTRAN is under the control of the IBJOB Mon- 
itor. See the IBJOB Processor publication, refer- 
enced above, for a description of the minimum ma- 
chine requirements for FORTRAN IV. 



MAJOR B.EVISION (May, 1963) 

This publication, Form C28-6274-1, is a major revision of the 
previous edition, Form C28-6274. It makes that publication, and 
the Techiiical Newsletter to that publication, N28-0051, obsolete. 



Copies of this and other IBM publications can be obtained through IBM Branch Offices. 

Address comments concerning the content of this publication to: 

IBM Corporation, Programming Systems Publications, Dept. D91, PO Box 390, Poughkeepsie, N.Y. 
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CHAPTER 1: GENERAL PROPERTIES OF A FORTRAN SOURCE PROGRAM 



WRITING THE SOURCE PROGRAM 

The statements of a FORTRAN source program are 
normally written on a standard FORTRAN coding 
sheet. An example of a FORTRAN source program 
is given in Figure 1-1; the purpose of this program 
is to determine the largest value attained by a set 
of numbers, Aj. 

1. Columns 1-5 of the first line of a statement 
may contain a statement number that is less 
than 32, 768, to identify the statement. 
Blanks and leading zeros are ignored in these 
columns. 

2. Column 6 of the first line of a statement 
must be a blank or a zero. 

3. Columns 7-72 contain the actual FORTRAN 
statement. Blanks are ignored except in an 
H field of a FORMAT or alphameric argu- 
ment or a DATA statement. 

4. A statement m^ay be continued over as many 
as nineteen continuation cards. Any card 
with a non-blank, non-zero column 6 is con- 
sidered a continuation card. 

5. Cards with a C in column 1 are not processed 
by FORTRAN, and columns 2-72 may be used 
for commients. 

6. Columns 73-80 are not processed by FOR- 
TRAN and may be used for identification. 



The order of execution of the source state- 
ments is governed by the normal source pro- 
gram statement sequencing given in Appendix A. 



PUNCHING THE SOURCE PROGRAM 

The FORTRAN statements, prepared as above, are 
punched into the standard FORTRAN card in Figure 
1-2 for input to the computer. 

THE FORTRAN STATEMENTS 

The FORTRAN statements may be divided into the 
following groups: 

1. The arithmetic statement specifies a numer- 
ical or logical calculation. 

2. The control statements govern the flow of 
control in the program. 
The input/output statements provide the nec- 



3. 



4. 



5. 



essary input/output routines and the input/ 
output format. 

The subprogram statements enable the pro- 
grammer to define and use subprograms. 

The specification statements provide infer - 
mation about the constants and variables used 
in the program and provide information about 
storage allocation. 
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FORTRAN STATEMENT 
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IDENTIFICATION 
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Figure, 1-2 



CHAPTER 2: CONSTANTS, VARIABLES, SUBSCRIPTS, AND EXPRESSIONS 



FORTRAN provides a means of expressing constants, 
variables, and a subscript notation for expressing 
one-, two-, and three-dimensional arrays of vari- 
ables. 

CONSTANTS 

Five types of constants are perm.itted in a FOR- 
TRAN source program: integer or fixed point, real 
or single-precision floating point, double-precision 
floating point, complex, and logical. 

Integer Constants 



Double- Precision Constants 



General Form 



An integer constant consists of 1-11 decimal digits written 
without a decimal point. 



Examples: 

3 

528 

8085 
An integer constant may be as large as 2^^-l, ex- 
cept when used for the value of a subscript or as an 
index of a DO, in which case the value of the integer 
is computed modulo 2 15, 



Real Constants 



General Form 


1. 


A real constant consists of 1-9 


significant decimal digits 




written with a decimal point. 






2. 


A real constant may be follow 


ed by a decimal exponent, 




which is written as the letter E followed 


by an integer 




constant. The field following 


the letter 


E must not be 




blank; it may be zero. 







Examples: 
21. 
.203 
8.0067 

5.0 E3 (means 5.0 x 10^, i.e. , 5000.) 
5.0 E-3 (means 5.0 x lO'^, i.e. , .005) 

1. The magnitude of a real constant must be be- 
tween the approximate limits of 10^° andl0~38^ 
or be zero. 

2. A real constant has precision to 8 digits. 



General Form 



1. A double -precision constant consists of 1-17 significant 
decimal digits written with a decimal point. 

2. If it is desired to specify a decimal exponent, or if the 
constant contains fewer than 10 digits, the letter D, 
followed by the exponent, must follow the number. The 
exponent is an integer constant. The field following the 
letter D must not be blank; it may be a zero. 



Examples: 

21.987538294 

21.9D0 

. 203D0 

5.0D3 (means 5.0x10^, i.e. ,5000.) 

5. OD-3 (means 5.0x10-3, i.e.,. 005) 

1. The magnitude of a double-precision constant 
must lie between the approximate limits of 
10-29 Qj^(^ 10^8, or be zero. Numbers be- 
tween lO"^^ and 10-38 j^ay i3g used, but only 
eight digits are significant in this range. 

2. Double-precision constants are floating point 
quantities that have precision to 16 digits. 

Complex Constants 



General Form 



A complex constant consists of an ordered pair of real constants 
separated by a comma and enclosed in parentheses. 



Examples: 

(3.2,1.86) is equal to 3.2+1.861. 

(2.1,0.0) is equal to 2.1 +<).0i. 

(5. 0E3, 2. 12) is equal to 5000. +2. 12i. 

1. The first real constant represents the real 
part of the complex number, and the second 
real constant represents the imaginary part 
of the complex number. 

2. The parentheses are required regardless of 
the context in which the complex constant 
appears. 

3. Each part of the complex constant may be 
preceded by a + or - sign. 

Logical Constants 



General Form 



A logical constant may take either of the following forms: 
.TRUE. 
. FALSE. 
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VARIABLES 



Form of Subscripts 



A variable is specified by its name and its type. 
There are five types of variables: integer, real, 
double-precision, complex, and logical. 

Variable Names 



General Form 



A variable name consists of 1-6 alphameric characters, the 
first of which must be alphabetic. 



Examples: 

L5 

JOBl 

BETATS 

COST 

K 
Subroutines are named in the same manner as 
variables (see "Naming Subroutines"). 



Variable Typ e Specification 

The type of a real or integer , variable or function 
name may be specified in one of two ways: implic- 
itly by name, or explicitly by a type statement (see 
the sections, "Type Statements" and "Naming Sub- 
routines"). All other variables must have their type 
specified by a type statement. 



General Form 


A subscript may 


take only 


one of the follow 


ing forms, 


where 


V represents any 


unsigned, 


non-subs. 


cripted 


integer variable, 1 


and c and c' represent any 

V 


unsignec 


integer 


constant: 






c 

v+c 












v-c 












c*v 












c*v+c' or c*v-c' 









Examples: 
IMAS 
J9 
K2 
N+3 

8*IQUAN 
5*L+7 
4*M-3 

7+2*K invalid 
9+ J invalid 

Subscripted Variables 



General Form 



A subscripted variable consists of a variable name followed by 
parentheses enclosing one, two, or three subscripts that are 
separated by commas. 



Implicit Type Assignment 

Implicit type assignment pertains only to integer 
and real, variable and function names, as follows: 

1. If the first character of the symbol is I, J, K, 
L, M, or N, it is an integer name; e.g. , 
MAX, JOB, IDIST, LESL. 

2. If the first character of the symbol is not I, 
J, K, L, M, or N, it is a real name; e.g. , 
ALPHA, BMAX, Q, WHIT. 



SUBSCRIPTS 

A variable may be made to represent any element of 
a one-, two-, or three-dimensional array of quan- 
tities by appending one, two, or three subscripts, 
respectively, to the variable name. The variable 
is then a subscripted variable. The subscripts are 
expressions of a special form whose value determines 
the member of the array to which reference is made. 



Examples: 
A(I) 
K(3) 

BETA (8*J+ 2, K-2,L) 
MAX (I, J, K) 

1. During execution, the subscript is evaluated 
so that the subscripted variable refers to a 
specific member of the array. 

2. Each variable that appears in subscripted 
form must have the size of the array speci- 
fied preceding the first appearance of the 
subscripted variable in an executable state- 
ment or DATA statement. This must be done 
by a DIMENSION or COMMON, with dimen- 
sions, statement. 

Arrangement of Arrays in Storage 

Arrays are stored in column order in increasing 
storage locations, with the first of their subscripts 
varying most rapidly and the last varying least rap- 
idly. 



Example: The 2-dimensional array Aj^, n is stored 
as follows, from lowest core storage location to 
highest: 

Al, 1, A2, 1, . . . , Am, 1. ^1, 2. A2, 2, • • • ' 
Ani,2» • • • . Ani,n 

EXPRESSIONS 

The FORTRAN language includes two kinds of ex- 
pressions: arithmetic and logical. 

Arithmetic Expressions 

An arithmetic expression consists of certain se- 
quences of constants, subscripted and nonsub- 
scripted variables, and arithmetic function refer- 
ences separated by arithmetic operation symbols, 
commas, and parentheses. 

The arithmetic operation symbols +, -, *, /, ** 
denote addition, subtraction, multiplication, divi- 
sion, and exponentiation, respectively. 

The following are the rules for constructing arith- 
metic expressions: 

1. Figures 2-1 and 2-2 indicate which constants, 
variables , and functions may be combined by the 
arithmetic operators to form arithmetic expressions. 
Figure 2-1 gives the valid combinations with respect 
to the arithmetic operators +, -, *, and/. Figure 
2-2 gives the valid combinations with respect to the 
arithmetic operator **. In these figures , Y indicates 
a valid combination and N indicates an invalid combi- 
nation. 



+,-,*,/ 


Real 


Integer 


Complex 


Double - 
Precision 


Logical 


Real 


Y 


N 


Y 


Y 


N 


Integer 


N 


Y 


N 


N 


N 


Complex 


Y 


N 


Y 


N 


N 


Etouble- 
Precision 


Y 


N 


N 


Y 


N 


Logical 


N 


N 


N 


N 


N 



Figure 2-1 



Exponent 



** 


Real 


Integer 


Complex 


Double - 
Precision 


Logical 


Real 


Y 


Y 


N 


Y 


N 


Integer 


N 


Y 


N 


N 


N 


Complex 


N 


Y 


N 


N 


N 


Double - 
Precision 


Y 


Y 


N 


Y 


N 


Logical 


N 


N 


N 


N 


N 



Base 



Figure 2-2 



2. A real constant, variable, or function name 
combined with a double-word quantity results in an 
expression with the type of the double-word quantity; 
e. g. , a real variable plus a complex variable forms 
a complex expression. 

3. Any expression may be enclosed in parenthe- 
ses. 

4. Expressions may be connected by the arithme- 
tic operation symbols to form other expressions, 
provided that: 

a. No two operators appear in sequence, and 

b. No operation symbol is assumed to be pres- 
ent. 

The expression A**B**C is not permitted; it must be 
written as either A**(B**C) or (A**B)**C, whichever 
is intended. 

5. Preceding an expression by a + or - sign does 
not affect the type of the expression. 

6. Hierarchy of Operations. Parentheses may be 
used in expressions to specify the order in which op- 
erations are to be computed. Where parentheses 
are omitted, the order is understood to be as fol- 
lows (fromi innermost operations to outermost opera- 
tions): 

Function Reference 
** Exponentiation 

* and / Multiplication and Division 

+ and - Addition and Subtraction 

For example, the expression A+B/C-D**E*F-G 
will be taken to mean A+ (B/C) - (dE*F)-G. 



a. 
b. 
c. 
d. 



Logical Expressions 

A logical expression consists of certain sequences of 
logical constants , logical variables , references to 
logical functions, and arithmetic expressions (ex- 
cept complex expressions) separated by logical op- 
eration symbols or relational operation symbols. A 
logical expression always has the value .TRUE, or 
. FALSE. . 

The logical operation symbols (where a and b are 
logical expressions) are: 



Symbol 
.NOT. a 



Definition 



This has the value . TRUE, only if 
a is . FALSE. ; it has the value 
. FALSE, only if a is .TRUE. 

a. AND.b This has the value .TRUE, only if 
a and b are both . TRUE. ; it has 
the value . FALSE, if either a or b 
is . FALSE. 

a. OR.b (Inclusive OR) This has the value 

. TRUE, if either a or b is .TRUE, 
it has the value . FALSE, only if 
both a and b are . FALSE. 
The logical operators NOT, AND, and OR must 
always be preceded and followed by a period. 
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The relational operation symbols are: 



Symbol 

7gt7~" 

,GE. 
.LT. 
.LE. 
• EQ. 

.NE. 



Definition 



Greater than 

Greater than or equal to 

Less than 

Less than or equal to 

Equal to 

Not equal to 
The relational operators must always be pre- 
ceded and followed by a period. 

The following are the rules for constructing logi- 
cal expressions: 

1. Figure 2-3 indicates which constants, vari- 
ables , and functions may be combined by the rela- 
tional operators to form a logical expression. In 
this figure, Y indicates a valid combination and N 
indicates an invalid combination. 

The logical expression will have the value . TRUE, 
if the condition expressed by the relational operator 
is met; otherwise, the logical expression will have 
the value .FALSE.. 

2. A logical expression may consist of a single 
logical constant, a logical variable, or a reference 
to a logical function. 



3. The logical operator .NOT. must be followed 
by a logical expression, and the logical operators 
.AND. and .OR. must be preceded and followed by 
logical expressions to form more complex logical 
expressions. 

4. Any logical expression may be enclosed in 
parentheses; however, the logical expression to 
which the . NOT. applies must be enclosed in paren- 
theses if it contains two or more quantities. 

5. Hierarchy of Operations. Parentheses may be 
used in logical expressions to specify the order in 
which operations are to be computed. Where paren- 
theses are omitted, the order is understood to be 

as follows (from innermost operation to outermost 
operation): 



Function Reference 

** 

* and / 

+ and - 

.LT., .LE., .EQ,,. 

.NOT. 

.AND. 

.OR. 



Exponentiation 
Multiplication and Division 
Addition and Subtraction 
NE., .GT.,.GE. 



.GT.,.GE.,.LT., 
.LE.,.EQ.,.NE. 


Real 


Integer 


Complex 


Double - 
Precision 


Logical 


Real 


Y 


N 


N 


Y 


N 


Integer 


N 


Y 


N 


N 


N 


Complex 


N 


N 


N 


N 


N 


Double - 
Precision 


Y 


N 


N 


Y 


N 


Logical 


N 


N 


N 


N 


N 



Figure 2-3 
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CHAPTER 3: THE ARITHMETIC STATEMENT 



The arithmetic statement defines a numerical or 
logical calculation. A FORTRAN arithmetic state- 
ment closely resembles a conventional arithmetic 
formula; however, the equal sign of the FORTRAN 
statement specifies replacement rather than equiva- 
lence. 



E and F are complex variables, and G, H, and P 
are logical variables. 



General Form 


a=b 










where: 










1. a is a real, integer, 


double 


-precision, 


complex, 


or 


logical, subscripted 


or non- 


-subscripted 


variable, 


and 


2. b is an expression. 











Examples: 

Ql =K 

A(I) = B(I) + ASIN(C(I)) 

V = .TRUE. 

E = C.GT.D.AND.F.LE.G 
Figure 3-1 indicates which type expressions may 
be equated to which type of variable in an arithmetic 
statement. In this figure, Y indicates a valid state- 
ment and N indicates an invalid statement. 







Right . 


ide of equ 


al sign 








\expression 
variable^y 


Real 


Integer 


Complex 


Double- 
Precision 


Logical 


Left 


Real 


Y 


Y 


N 


Y 


N 


side 


Integer 


Y 


Y 


N 


Y 


N 


of 


Complex 


N 


N 


Y 


N 


N 


equal 
sign 


Double - 
Precision 


Y 


Y 


N 


Y 


N 




Logical 


N 


N 


N 


N 


Y 



Figure 3-1 

In the following examples of arithmetic state- 
ments, I is an integer variable, A and B are real 
variables, C and D are double -precis ion variables, 



A = B 
I = B 



A = I 



1 = 1+ 1 

A= 3*B 



A = B*C 



E=F*(3.7,2.0) 



G = .TRUE. 
H= .NOT.G 



H = I.GE.A 



G = H.OR..N0T.P 



G= 3..GT.B 



Replace A by the current value of B. 

Truncate B to an integer, convert it to 

an integer constant, and store it in I. 

Convert I to a real variable and store 

it in A. 

Add 1 to I and store it in I. 

Not permitted. The expression is mixed 

for multiplication, i.e., it contains 

both a real variable and an integer 

constant. 

Multiply B by C using double-precision 

arithmetic, and store the most significant 

part of the result as a real number in A. 

Multiply F by 3.7 + 2.0i using complex 

arithmetic, and store the result in E as 

a complex number. 

Store the logical constant .TRUE, in G. 

If G is .TRUE., store the value .FALSE., 

in H; if G is . FALSE. , store the value 

• TRUE, in H. 

Not permitted. An integer and a real 

variable may not be joined by a relational 

operator. 



where: 
-^ implies 
• NOT., asd 
V implies 
.OR. 



G is .TRUE, if 3. is greater than B; G is 
.FALSE, otherwise. 



H 


P 


— P 


Hv — P 


T 


T 


F 


T 


T 


F 


T 


T 


F 


T 


F 


F 


F 


F 


T 


T 



The last two examples above illustrate the follow-: 
ing rules: 

1. Two logical operators may appear in se- 
quence, not separated by a comma or paren- 
theses, only if the second logical operator 

is .NOT. . 

2. Two decimal points may appear in succession 
as in item 1 or when one belongs to a con- 
stant and the other to a relational operator. 



The Arithmetic Statement 11 



CHAPTER 4: THE CONTROL STATEMENTS 



The control statements enable the programmer to 
control and terminate the flow of his program. 

The Unconditional GO TO Statement 



The ASSIGN Statement 



General Form 



GO TOn 

where: 

n is a statement number. 



Example: 

GO TO 25 
This statement causes control to be transferred 
to the statement numbered n. 

The Computed GO TO Statement 



General Form 



GOTO(ni, n2, ...,nm),i 
where: 

1. nj, n2, ... ,njjj are statement numbers, and 

2. i is a non-subscripted integer variable. 



Example: 

GO TO (30, 45, 50, 9), K 
This statement causes control to be transferred 
to the statement numbered ni, n2, . . . , n-^ depending 
on whether the value of i is 1, 2, 3, . , , m, respec- 
tively, at the time of execution. Thus, in the ex- 
ample, if K is 3 at the time of execution, a transfer 
to the third statement in the list, i.e. , statement 
50, will occur. 

The Assigned GO TO Statement 



General Form 



GO TOn, (ni, ng, ...,nm) 
where: 

1. n is a non-subscripted integer variable appearing in a 
previously executed ASSIGN statement, and 

2. nj, n2, i, . . , nm are statement numbers. 



Example : 

GO TO J, (17,12,19) 
This statement causes control to be transferred to 
the statement number last assigned to n by an AS- 
SIGN statement; n^, n2, . . . .nj^ is a list of the m 
values that n may assume. 



General Form 



ASSIGN i TO n 

where: 

1 . i is a statement number, and 

2. n is a non-subscripted integer variable that appears in 
an assigned GO TO statement. 



Examples : 

ASSIGN 12 TO K 

ASSIGN 37 TO JA 
This statement causes a subsequent GO TO n, (ni, 
n2, . . . ,njn) to transfer control to the statement num- 
bered i, where i is one of the statement numbers in- 
cluded in the series ni, n2, ... jnm. 

The Arithmetic IF Statement 



General Form 



IF (a) ni,n2,n3 
where: 

1. a is an arithmetic expression (not complex), and 

2. nj,n2,n3 are statement numbers. 



Examples : 

IF (A(J,K) - B) 10,4,30 

IF (D*E+BRN) 9, 9, 15 
This statement causes control to be transferred to 
the statement numbered n^, n2, or n^ if the value of 
a is less than, equal to, or greater than zero, re- 
spectively. 

The Logical IF Statement 



General Form 



IF (t) s 

where: 

1. t is a logical expression, and 

2. s is any executable statement except DO or another 
logical IF. 



Examples: 

IF (A. AND. B) F = SIN (R) 
IF (16.GT.L) GO TO 24 
IF (D.OR.X.LE.Y) GO TO (18, 20), I 
IF (Q) CALL SUB 
1. If the logical expression t is true, statement 
s is executed. Control is then transferred to 
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the next sequential statement unless s is an 
arithmetic IF-or GO TO-type statement, in 
which case, control is transferred as indicated. 
If t is false, conti'ol is transferred to the next 
sequential statement. 

If t is true and s Is a CALL statemtmL, upon 
return from the subprogram control is trans- 
ferred to the next sequential statement. 



The DO Statement 



General Form 



DO n i = m^, m^} mj 

where: 

1. n is a statement number, 

2. i is a non-subscripted inte'ger variable, and 

3. mijniojnig are each either an unsigned integer constant or anon 
subscripted integer variable; if m^^ is not stated, it istakento be 1. 



Examples; 

DO 30 I - 1, M, 2 
DO 24 I = 1, 10 

The DO statement is a command to execute re- 
peatedlj' the statements that follow, up to and in- 
cluding the statement numbered n. The statements 
in the range of the DO are executed repeatedly with 
i equal to' mi) then i equal to m]_ + ms, then i equal 
to mi + 2mo, etc. , until i is equal to the higb'^st 
value in this sequence that does not exceed rw^- The 
state;>^ent8 in the range of the BO will be cxccuiod 
at least once. 

1.. The range of a DO is that set of statements 
that will be executed repeatedly; i.e. , it is tlie se- 
quence of consecutive statements immediately follow 
Ing the DO statement, up to and including the state- 
ment numbered n. After the last cxecutioii of ir.t; 
range, the DO is said to be satisfied. 

2. The index of a DO is the integer variable i. 
Tliroughout the range of the DO, the index is avail- 
able for computation, either as an ordinary integer 
variable or as the variable of a subscript, uoo.i 
exiting from a DO by satisfying the DO, the index 1 
must be redefined before it is used in computation. 
Upon exiting from a DO by transferring out of the 
range of the DO, the index i is available for com- 
putation and is equal to the last value it attained. • 

3. DOs within DOs . Among the statements in the 
range of a DO may be other DO statements; such a 
configuration is called a nest of DOs. If the nuigc 
of a DO includes another DO, then all of the state- 
ments in the range of the latter must also be in the 
range of the former. 

4. Transfer of Control and DOs . Control may not 
be transferred into the rtrnge of a DO from outside 
its rimge. Thus, in the configuration below 1,2, and 
3 are permitted transfers, but 4,5, and 6 are not. 
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DO 



4 X 




5. R estrictions on Statements in the Rang e of a 
DO 

a. Any statement that redefines the index or 
any of the indexing parameters (m's) is 
not permitted in the range of a DO. 

b. The range of a DO cannot end with an arith- 
metic IF or GO TO-type statement or with 
a non-executable statement. The rAnge of a 
DO may end with a logical IF, in which case, 
control is handled as follows: if the logical 
expression t is false, the DO is reiterated; 
if the logical expression t is true, -statement 
s is executed and then the DO is reiterated. 
However, if t is true and s is an arithmetic 
IF or transfer type statement, control is 
transferred as indicated. 

6. vViion a reference to a subprogram is executed 
in the range of a DO, care must be taken that the 
called subprogram does not alter the DO index or 
the indexing parameters. 

The CONTINUE Statement 



General Form 



CONTINUE 



CONTINUE is a dummy statement that givos rise 
to no instructions in the object program. It is most 
frequently used as the last statement in the range of 
a DO to provide a transfer address for IF and GO TO 
statements that are intended to begin another repe- 
tition of the DO range. 
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The PAUSE SluLomcnL 



The END statement must be the physically 
last statement of the program. 



General Form 



The STOP Statement 



PAUSE or PAUSE n 

where n is an unsigned octal integer constant of 1-5 digits. 



Examples : 

PAUSE 

PAUSE 77777 
The machine will halt with the octal number n 
ll>u uddi L'Syfreiiil ul LliL iii i t i Diit>g^i» > .3B i< jiig;r ' j!rti!jii>r If n is noi 
specified, it is understood to be zero. Depressing 
the START key causes the program to resume exe- 
cution of the object program with the next executable 
FORTRAN statement. 

The END Statement 



Ac 



General Form 



END 



1. The END statement terminates compilation 
of a program. 



General Form 



STOP 



The STOP statement terminates the execution of 
any program by returning control to the Monitor. 
Execution of a program may also be terminated by 
a CALL to the EXIT and DUMP subroutines (see 
the section, "EXIT, DUMP, and PDUMP"). 

The RETURN Statement 



! General Form 



RETURN 



The normal exit from any subprogram is the 
P.ETURN statement, which returns control to the 
calling program. The RETURN statement is the 
logical end of the program; there may be any num- 
ber of RETURN statements in the program. 
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CHAPTER 5: INPUT/OUTPUT STATEMENTS 



The FORTRAN statements that specify transmission 
of information to or from I/O devices may be grouped 
as follows: 

General I/O Statements: The statements READ and 
WRITE cause the transmission of a specified list of 
quantities between core storage and an input/output 
device. The statements PUNCH and PRINT cause 
information to be transmitted from core storage to 
the card punch and on-line printer, respectively. 

Manipulative l/O Statements: Statements END FILE , 
REWIND, and BACKSPACE manipulate l/O devices. 

Format Statement: The non-executable statement 
FORMAT specifies the arrangement of the informa- 
tion in the external input/output medium specified 
by the general I/O statements. 

LIST SPECIFICATIONS 

The general I/O statements must include an ordered 
list of the quantities to be transmitted. The order of 
the I/O list must be the same as the order in which 
the information exists in the input/output medium. 

The following notes on the formation and meaning 
of an I/O list are most clearly understood by con- 
sidering the following I/O list: 

A,B(3), (C(I), D(I,K), 1=1, 10), 
((E(I,J), 1=1, 10,2), F(J,3), J=1,K) 
which implies that the information in the external 
I/O medium is arranged as follows: 

A,B(3),C(1),D(1,K), C(2), D(2,K),..., 
C(10), D(10,K), E(l,l), E(3,l),..., 
E(9,l), F(l,3), E(l,2), E(3,2), .... 
E(9,2), F(2,3) F(K,3) 

1 . An I/O list is a string of list items separated 
by commas. A list item may be: 

a. A subscripted or non-subscripted variable; 
or 

b. An implied DO. 

An I/O list reads from left to right, with repetition 
of variables enclosed in parentheses. 

2. A constant may appear in an I/O list only as a 
subscript or as an indexing parameter. 

3. The execution of an I/O list is exactly that of 
a DO-loop, as though each left parenthesis (except 
subscripting parentheses) were a DO, with indexing 
given immediately before the matching right paren- 
thesis , and with the DO range extending up to that 
indexing information. The order of the l/O list 
above may be considered equivalent to the following 
"program": 



(C(I),D(I,K),I=1,10) 



((E(I,J),I=1,10,2), 
F(J,3),J=1,K) 



A 

B(3) 

DO 5 I = 1,10| 
C(I) I 

5 D(I,K) ^ 

DO 9 J=1,K 
DO 8 1=1,10,21 

8 E(I,J) 

9 F(J, 3) 

4. An implied DO is best defined by an example. 
In the I/O list above, the list item (C(I),D(I,K), 
1=1, 10) is an implied DO; it is evaluated as in the 
above program. 

The range of an implied DO must be clearly de- 
fined by parentheses. 

5. For a list of the form K, A(K), or K, (A(I), 
1=1, K), where the definition of an index or an index- 
ing parameter appears earlier in the list of an input 
statement than its use, the indexing will be carried 
out with the newly read in value. 

6. Any number of quantities may appear in a sin- 
gle list. However, each quantity must have the cor- 
rect format as specified in a corresponding FOR- 
MAT statement. Essentially, it is the list which 
controls the quantity of data read. If more quan- 
tities are to be transmitted than are in the list, only 
the number of quantities specified in the list are 
transmitted, and remaining quantities are ignored. 
Conversely, if a list contains more quantities than 
are given on one BCD input record, more records 
are read; if a list contains more quantities than are 
given in one binary record, reading is terminated 

as an object program error and control is trans- 
ferred to FXEM. 

INPUT/OUTPUT OF ENTIRE ARRAYS 

When input/output of an entire array is desired, an 
abbreviated notation may be used in the list of the 
input/output statement; only the name of the array 
need be given, and the indexing information may be 
omitted. 

1. If A has previously been listed in a DIMEN- 
SION or COMMON (with dimensions) state- 
ment, the statement 
READ (5, 10) A 
is sufficient to read in all of the elements of 
the array A (see the section, "Input"). 
2. The elements read in by this notation are 
stored in accordance with the description of 
the arrangement of arrays in storage (see 
the section, "Arrangement of Arrays in Stor- 
age"). 
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3. If A has not been previously dimensioned, 
only the first element will be read in. 

FORMAT 

The BCD input/output statements require, in addition 
to a list of quantities to be transmitted, reference to 
a FORMAT statement which describes the type of 
conversion to be performed between the internal ma- 
chine language and the external notation for each 
quantity in the list. 



General Form 



FORMAT (SpS2,...,Sn/S'i,S'2,.. 
field, Sj, is a format specification. 



, S'n/ . . . ) where each 



Example : 

FORMAT (12/ (E12.4,F10.2)) 

1. FORMAT statements are not executed; they 
may be placed anywhere in the source program. 
Each FORMAT statement must be given a statement 
number. 

2. The FORMAT statement indicates, among 
other things, the maximum size of each record to 
be transmitted. In this connection, it must be re- 
membered that the FORMAT statement is used in 
conjunction with the list of some particular input/ 
output statement, except when a FORMAT statement 
consists entirely of alphameric fields. In all other 
cases, control in the object program switches back 
and forth between the list (which specifies whether 
data remains to be transmitted) and the FORMAT 
statement (which gives the specifications for trans- 
mission of that data). 

3. Slashes are used to specify unit records, 
which must be one of the following: 

a. A tape record with a maximum length 
corresponding to the printed line of the 
off-line printer. 

b. A punched card, to be read on-line, with 
a maximum of 72 characters; a punched 
card, to be read off-line, with a maxi- 
mum of 80 characters. 

c. A line to be printed on-line, with a maxi- 
mum of 120 characters. 

Thus, FORMAT (3F9.2,2F10.4/8E14.5) would speci- 
fy records in which the first, third, fifth, etc. , have 
the format (3F9.2, 2F10.4), and the second, fourth, 
sixth, etc. , have the format (8E14.5). 

4. During input/output of data, the object pro- 
gram scans the FORMAT statement to which the 
relevant input/output statement refers. When a 
specification for a numerical field is found and list 
items remain to be transmitted, input/output takes 
place according to the specification, and scanning 
of the FORMAT statement resumes. If no items 



remain, transmission ceases and execution of that 
particular input/output statement is terminated. 
Thus, a decimal input/output operation will be 
brought to an end when there are no items remaining 
in the list. 

Numeric Fields 

Five types of conversion are available for numeric 
data: 



Internal 


Conversion Code 


External 


Floating point 

( double -precision) 


D 


Real with D exponent 


Floating point 


E 


Real with E exponent 


Floating point 


F 


Real without exponent 


Integer 


I 


Integer 


Integer 


O 


Octal Integer 



These types of conversion are specified in the 
forms Dw.d, Ew.d, Fw.d, Iw, Ow, where: 

1. D, E, F, I, and O represent the type of con- 
version; 

2. w is an unsigned integer constant that repre- 
sents the field width for converted data; this 
field width may be greater than required in 
order to provide spacing between numbers; 

3. d is an unsigned integer or zero that repre- 
sents the number of positions of the field that 
appear to the right of the decimal point. 

For example, the statement FORMAT (I2,E12.4, 
OB, F10.4, D25.16) might cause the following line to 
be printed: 

I2Ei2.4 08 F10.4 D25. 16 

27lb-0. 9321Eb02(57734276bbb-0. 0076bb-0. 7878977909500672Db03 

where b indicates a blank space. 

Notes on D-, E-, F-, I-, and 0-Conversion 

1. Specifications for successive fields are sep- 
arated by commas . 

2. No format specification that provides for more 
characters than permitted for a relevant input/out- 
put record should be given. Thus, a format for a 
BCD record to be printed off-line should not provide 
for more characters (including blanks) than the capa- 
bilities of the printer. 

3. Information to be transmitted with 0-conversion 
may be given any type of variable name; information 
to be transmitted with E- and F-conversion must 
have real names; information to be transmitted with 
I-conversion must have integer names; information 
to be transmitted with D-conversion must have 
double-precision names. 



16 



4. The field width w, for D-, E-, and I^ -conver- 
sion, must include a space for the decimal point and 

a space for the sign. Thus, for D- and E-co!v,Trsion, 
\v>di7, and for F-convcrsion, Wi^d+S. 

5. The exponent, which may be used with D- and 
E-conversion, is t!ic power of 10 to which the ivum- 
ber must be raised to obtain its true value. 1".- ex- 
ponent is written with an E (for E-conversiosi) or 1) 
(for D-conversion) followed by a minus sign if the 
exponent is negative, or a plus sign or a bianlv if the 
exponent is positive, and then followed by two num- 
bers which arc the exponent. For example, the num- 
ber .002 is equivalent to the number ,2E-02. 

6. For input under D-conversion, up to 17 deci- 
mal digits arc converted and the result is stored so 
that the most significant part and the least signifi- 
cant part are in adjacent core storage locations. 

For output under D-conversion, the two core stor- 
age words representing the double-precision quan- 
tity arc considered one piece of data and converted 
as such.' 

7. If a number converted by I-conversion re- 
quires more spaces than are allowed by tlie field 
width w, the excess on the high-order side is lost. 
If the number requires fewer than w spaces , the 
leftmost spaces are filled with blanks. If the number 
is negative , the space preceding the leftmost digit 
will contain a minus sign if sufficient spaces have 
been reserved. 

Complex Nu mber Fields 

Since a complex quantity consists of two separaic 
and independent real numbers, a complex number 
is transmitted by two successive real number speci- 
fications or by one real number specification that is 
repeated. 

Sec Figure Tr-l for an example of a FORMAT 
statement to transmit an array consisting of six 
complex numbers. 



FORMAT (2E10. 2, E8.3, 1PE9.4, E10.2, F8. 4, 3(EiO, 2, fa . 2)) 



Figure 5-1 

Alphameric Fields 

FORTRAN provides two ways by which alphameric 
information may be transmitted; both specifications 
result in storing the alphameric information inter- 
nally in BCD, 

1. The specification Aw causes w characters to 
be read into, or written from, a varial)le or 
array name. 

2. The specification nil introduces alphameric 
information into a FORMAT statement. 



The basic difference between A- and li-conversion 
is that information handled by A-conversion is given 
a variable name or array name and hence can be re- 
ferred to by means of 'this name for processing and 
modification, whereas, information handled by H- 
conversion is not given a name and may not be re- 
ferred to or manipulated in storage in any way, 

A-Conversion 

The variable name to be converted by A-conversion 
must conform to the normal rules for naming FOR- 
TRAN variables; it may be any type. 

1. On input, nAw will be interpreted to mean 
that the next n successive fields of w charac- 
ters each are to be stored as BCD informa- 
tion. If w is greater than 6, only the 6 right- 
most characters will be significant. If w is 
less than 6, the characters will be left-ad- 
justed, and the word filled out with blanks. 

2, On output, nAw will be interpreted to moan 
that the next n successive fields of w char- 
acters each are to be the result of transmis- 
sion from storage without conversion. If w 
exceeds 6, only 6 characters of output will 
be transmitted, preceded by w-6 blanks. If 
w is less than 6, the w leftmost characters 
of the word will be transmitted, 

H-Conversion 

The specification nil is followed in the FORMAT 
statement by n alphameric characters. For example, 

31H THIS IS ALPHAMERIC INFORMATION 
Note that blanks are considered alphameric charac- 
ters and must be included as part of the count n. The 
effect of nH depends on whether it is used with input 
or output, 

1, On input, n characters are extracted from 
the input record and replace the n characters 
included with the source program FORMAT 
specification. 

2, On output, the n characters following the 
specification, or the characters that replaced- 
them, are written as part of the output rec- 
ord. 

See Figure 5-2 for an example of A- and H-con- 
version in a FORMAT statement. 

The statement FORMAT (4HbXY-, F8.3, A8) might 
produce the following lines , where b indicates a 
blank character: 



XY 


= b-93 


210bbbbbbbb 


XY 


= 9999 


999bbOVFLOW 


XY 


= bb28 


768bbbbbbbb 



Figure 5-2 
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Figui-c 5-2 assumes that there arc steps in the 
source program that read the data OVFLOW, store 
this data in the word to be printed in the for!!] at A^ 
when overflow occurs, and store six blanks in the 
word when overflow does not occur. 

Logical Fields 

Logical variables may be read or written by means 
of the specification Lw, where L rejoresents the logi- 
cal type of conversion and w is an integer constant 
that represents the data field width. 

1, On input, a value of either ti-ue or false will 
be stored if the first non-blank character in 
the field of w characters is a T or an F, re- 
spectively. If all'the w characters are blank, 
a value of false will be stored. 

2. On output, a value of true or false in storage 
will cause w minus 1 blanks, followed bj^ a T 
or an F, respectively, to be written out. 

Blank Fields - X-Conversion 



The specification nX introduces n blank characters 
into an input/output record where < n < 132. 

1. On input, nX causes n characters in the input 
record to be skipped, regardless of what they 
actually are. 

2. On output, nX causes n blanks to be intro- 
duced into the output record. 

Repetition of Field Format 

It may be desired to print or read n successive fields 
in the same format within one record. This may be 
specified by giving n, an unsigned integer, before 
D,E,F,I,L,0, or A. Thus, the field specification 
aE12.4 is the same as writing E12,4, E12.4, E12.4. 

Repetition of Groups 

A limited parenthetical expression is permitted in 
order to enable repetition of data fields according 
to certain format specifications within a longer FOR- 
MAT statement. Thus, FORMAT (2(F10.6,E10,2),I4) 
i;3 equivalent to FORMAT (FlO. 6, ElO. 2, FIO. 6, 
E10.2, 14). Two levels of parentheses, in addition 
to the parentheses required by the FORMAT state- 
ment, are permitted. The second level of paren- 
theses facilitates the transmission of complex quan- 
tities. 

Scale Factors 



To permit more general use of D-, E-, and F-con- 
version, a scale factor followed by the letter P may 



prccc(!e the specification. The magnitude of the 

scale factor must be between -8 and ^•8, inclusive. 

Tho scale factor is defined for input as follows: 

2()-scalc factor x external quantity == internal 

quantity 

The scale factor is defined for output as follows; 

external quantity^intemal quantity 
^ ] Q scale factor 

For input, scale factors have -effect only on F- 
conversion. For example, if input data is in the 
form xx.xxxx mid it is desired to use it internally in 
the form .xxxxxx, then the FORMAT specification 
to effect this change is 2PF7.4. For output, scale 
factors may be used with D-, E-, and F-conversion. 
For example, the statement FORMAT (I2,3F11.3) 
might give the following printed line: 

27bbbb-93. 209bbbbb-0. OOSbbbbbbO. 554 
But the statement FORMAT (12, 1P3F11. 3) used with 
the same data would give the following line: 

27bbb-932. 094bbbbb-0. 076bbbbbb5. 536 
Whereas, the statement' FORMAT (12, -1P3F11.3) 
wouLtl give the following line: 

v7bbbbb-9. 321bbbbb-0. OOlbbbbbbO. 055 
A positive scale factor used for outputwith D- and 
E-convcrsion increases the number and decreases the 
exponent. Thus, FORMAT (12, 1P3E12. 4) would 
produce with the same data: 

27b-9.3209Eb01b-7.5804E-03bb5.5536E-01 
The scale factor is assumed to be zero if no other 
value has been given. However, once a value has 
been given, it will hold for all D-, E-, and F-con- 
versions following the scale factor within the same 
FORMAT statement. This applies to both single- 
record formats and multiple -record formats (see 
below). Once the scale factor has been given, a 
subsequent scale factor of zero in the same FOR- 
MAT statement must be specified by OP. Scale fac- 
tors have no effect on I- and O-conversion. 

Multiple-Record Formats 

To deal with a block of more than one line of print, 
a FORMAT specification may have several different 
one-line formats separated by a slash,/, to indicate 
the iKiginning of a new blank line. Thus, FORMAT 
(3Fy.2,2F10.4/8E14.5) would specify a multi-line 
block of print in which lines 1, 3, 5, ... have format 
(3F9.2,2F10.4), and lines 2, 4, 6, ... have format 
(8E14.5). 

If a m.ultiple-line format is desired in which the 
first two lines are to be printed according to a spe- 
cial format and all remaining lines according to an- 
other format, the last line-specification should be 
enclosed in a second pair of parentheses; e.g. , 
FORMAT (I2,3E12.4/2F10.3,3F9.4/ (10F12.4)). 
If data items remain to be transmitted after the for- 
mat specification has been completely "used," the 
format repeats from the last previous parenthesis 
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which is a zero or a first level parenthesis, 
example, consider the FORMAT statement: 



For 



FORMAT (3 El 0.3, (12, 2 (F12.4)), D28.17) 
1 2 21 

The parentheses labeled are level parentheses; 
those labeled 1 are first level parentheses; and, 
those labeled 2 are second level parentheses. If 
more items in the list are to be transmitted after 
the format statement has been completely used, the 
FORMAT repeats from the last first level left paren- 
thesis; i.e. , the parenthesis preceding 12. 

As these examples show, both the slash and the 
right parenthesis of the FORMAT statement indicate 
a termination of a record. 

Blank lines may be introduced into a multi-line 
FORMAT statement by listing consecutive slashes. 
When n+1 consecutive slashes appear at the end of 
the FORMAT, they are treated as follows: for input, 
n+1 records are skipped; for output, n blank lines 
are written. When n+1 consecutive slashes appear 
in the middle of the FORMAT, n records will be 
skipped for both input and output. 

Carriage Control 

The WRITE (i,n) list statement prepares a BCD tape 
that can be used to obtain off-line printed output. 
The PRINT n, list statement prints on-line during 
execution. The off-line printer may be set manually 
to operate in one of three modes: single space, 
double space, and program control. Under program 
control, the first character of each BCD record, 
the control character, controls spacing of the printer; 
this first character is not printed. The control 
characters and their effects, both on-line and off- 
line, are: 



DIMENSION FMT (12) 
FORMAT (12A6) 
READ (5, 1) (FMT(I), 1=1, 12) 
READ(5,FMT) A,B,(C(I), 1=1,5) 



Character 


Effect 


Off-line 


On-line 


Blank 



1 


Single space 
Double space 
Eject 


before printing 
before printing 
before printing 


after printing 
after printing 
before printing 
followed by a 
single space 
after printing 



Program control is usually obtained by beginning 
a FORMAT specification, for a BCD record, with 
IH followed by the desired control character. 

FORMAT Statements Read In at Object Time 

FORTRAN accepts a variable FORMAT address. 
This permits specifying a FORMAT for an I/O list 
at object time. 



Figure 5-3 

In Figure 5-3, A, B, and the array C are converted 
and stored according to the FORMAT specifications 
read into the array FMT at object time. 

1. The name of the variable FORMAT specifica- 
tion must appear in a DIMENSION statement, 
even if the array size is only 1. 

2. The format read in at object time must take 
the same form as a source program FORMAT 
statement, except that the word FORMAT is 
omitted, i.e., the variable format begins 
with a left parenthesis. 

Data Input To the Object Program 

Data input to the object program is punched into 
cards according to the following specifications: 

1. The data must correspond in order, type, 
and field with the field specifications in the 
FORMAT statement. Punching begins in 
card column 1. 

2. Plus signs may be omitted or indicated by a 
+. Minus signs are indicated by an 11 -punch. 

3. Blanks in numeric fields are regarded as 
zeros. 

4. Numbers for E- and F-conversion may con- 
tain any number of digits, but only the high- 
order 8 digits of precision will be retained. 
For D-conversion, the high-order 16 digits 
of precision will be retained. In both cases, 
the number is rounded to 8 or 16 digits of 
accuracy, as applicable. 

To permit economy in punching, certain relaxations 
in input data format are permitted. 

1. Numbers for D- and E-conversion need not 
have four columns devoted to the exponent 
field. The start of the exponent field must 
be marked by a D or E or, if that is omitted, 
by a + or - (not a blank). Thus, E2, E+2, 
+2, +02, and D+02 are all permissible ex- 
ponent fields . 

2. Numbers for D-, E-, and F-conversion need 
not have their decimal point punched; the for- 
mat specification will supply it. For exam- 
ple, the number -09321+2 with the specifica- 
tion E12.4 will be treated as though the deci- 
mal point had been punched between the and 
the 9. If the decimal point is punched in the 
card, its position overrides the position in- 
dicated in the FORMAT specification. 
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THE GENERAL INPUT/OUTPUT STATEMENTS 

Input 

The READ statement designates input. The follow- 
ing table gives the forms of the READ statement, 
where i, an unsigned integer constant or an integer 
variable, is a reference to an input device, and n 
is a FORMAT statement number. 



Output 

The statements PRINT and PUNCH are used to des- 
ignate the printing and punching of data on-line. All 
other output is designated by the statement WRITE. 
The following table gives the formats of the output 
statements, where i, an integer constant or an in- 
teger variable , is a reference to an output device , 
and n is a FORMAT statement number. 



Type of Input 


General Form 


Cards on-line 
BCD record 
Binary record 


READ n, list 
READ (i,n) list 
READ (i) list 



Examples : 

READ 10, (A(I),I=1,5) 
READ (5,10) A, B, (D(J), J=l, 10) 
READ (N, 10) K,DC(J) 
READ (3) (A(J),J=1,10) 
READ (N) (A(J), J=l,10) 

1. The READ n, list statement causes cards to 
be read from the card reader. 

2. The READ (i,n) list statement causes BCD 
information to be read from symbolic input 
device i (except the card reader). 

3. The READ (i) list statement causes binary 
information to be read from symbolic input 
device i (except the card reader). 

4. Under the first two forms of the READ state- 
ment, successive records are read until the 
entire I/O list has been satisfied: i.e. , all 
data items have been read, converted, and 
stored in the locations specified by the I/O 
list. 

Under the form READ (i) list, a record is read 
completely only if the list specifies as many words 
as the record contains. Binary records to be read 
in by a FORTRAN program should be written by a 
FORTRAN program or should be in the proper bi- 
nary record format as follows: 

Consider a logical record as being any sequence 
of binary words to be read by any one input state- 
ment. This logical record must be divided into 
physical records, each of which is a maximum of 
256io words long. Of course, if a logical record 
consists of fewer than 256io words, it will com- 
prise only one physical record. The first word of 
each physical record is a "signal" word that is not 
part of the list. This word contains zero for all but 
the last physical record of a logical record. The 
first word of the last physical record contains a 
number designating the number of physical records 
in this logical record. 



Type of Output 


General Form 


Cards on-line 
Print on-line 
BCD tape record 
Binary tape record 


PUNCH n, list 
PRINT n, list 
WRITE (i,n) list 
WRITE (i) list 



Examples : 

PUNCH 20, (A(J), J=l,6) 
PRINT 2, (A(J), J=l,6) 
WRITE (6,10) A,B, (C(J),J=1,10) 
WRITE (N, 11) K, D(J) 
WRITE (2) (A(J), J=l, 10) 
WRITE (M) A, B, C 

1. The PUNCH statement causes alphameric 
cards to be punched on-line. 

2. The PRINT statement causes data to be out- 
put on the on-line printer. 

3. The WRITE (i,n) list statement causes BCD 
information to be written on symbolic output 
device i. 

4. The WRITE (i) list statement causes binary 
information to be written on symbolic output 
device i. 

5. In the output statements in items 1, 2, and 3, 
successive records are written in accordance 
with the FORMAT statement until the list has 
been satisfied. In the output statement in 
item 4, one logical record is written consist- 
ing of all the words specified in the list. 

THE MANIPULATIVE INPUT/OUTPUT STATE- 
MENTS 

The statements END FILE, REWIND, and BACK- 
SPACE manipulate I/O devices. In the following 
table, i, an unsigned integer constant or integer 
variable, is a reference to an input/output device. 



General Form 



END FILE i 
REWIND i 
BACKSPACE i 
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2. 



3. 



Examples : 

END FILE 3 

END FILE N 

REWIND 3 

REWIND N 

BACKSPACE 3 
1. The END FILE i statement causes an end- 
of-file mark to be written on symbolic tape i. 
The REWIND i statement causes symbolic 
tape unit i to be rewound. 
A request to write an end of file or rewind 
system files SYSINl, SYSOUl, and SYSPPl, 
corresponding in the standard FORTRAN I/O 
Library to symbolic units 5,6, and 7, will 
be ignored. 
4. The BACKSPACE i statement causes tape i 
to be backspaced one physical record if i re- 
fers to an I/O device in the BCD mode , or it 
causes tape i to be backspaced one logical 
record if i refers to an I/O device in the bi- 
nary mode. A request to backspace SYSOUl 
corresponding in the standard FORTRAN I/O 
Library to symbolic unit 6 , will be ignored. 

SYMBOLIC INPUT/OUTPUT UNIT DESIGNATION 

Input/output devices are always referred to sym- 
bolically in FORTRAN input/output statements. 



Object program input/output operates 
through the Minimum IOCS Buffering package. 
The correspondence between the symbolic 
unit reference and the actual physical unit is 
established in the initialization of IOCS. 
The standard FORTRAN I/O configuration 
allows for symbolic tape units 1 through 8. 
The normal unit designation for BCD input 
statements is 5; for BCD output statements 
it is 6. 



FORTRAN Tape Units 


Mode 


Function 


1 


Binary 


Input or Output 


2 


Binary 


Input or Output 


3 


Binary 


Input or Output 


4 


Binary 


Input or Output 


5 


BCD 


Input 


6 


BCD 


Output 


7 


Binary 


Output 


8 


BCD 


Input or Output 



The symbolic unit references may be changed by 
each installation in accordance with its own needs. 
See the reference manual IBM 7090/7094 Opera ting 
Systems: Basic Monitor (IBSYS), Form C28-6248. 
and the reference manual IBM 7090/7094 Prog ram- 
ming Systems, IRJOB Processor, Form C28-6275. 
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CHAPTER 6: SUBROUTINES, FUNCTIONS, AND SUBPROGRAM STATEMENTS 



There are four classes of subroutines in FORTRAN: 
Arithmetic Statement Functions, Built-in Functions, 
FUNCTION subprograms, and SUBROUTINE sub- 
programs. The major differences among the four 
classes of subroutines are as follows: 

1. The first three classes may be grouped as 
functions; they differ from the SUBROUTINE 
subprogram in the following respects: 

a. The functions are always single-valued 
(that is, they return a single result); 
the SUBROUTINE subprogram may 
return more than one value. 

b. A function is referenced by an arithme- 
tic expression containing its name; a 
SUBROUTINE subprogram is refer- 
enced by a CALL statement. 

2. The built-in function is an open subroutine; 
i.e. , a subroutine that is incorporated into 
the object program each time it is referred 
to in the source program. The three other 
FORTRAN subroutines are closed; i.e. , they 
appear only once in the object program. 

NAMING SUBROUTINES 

All four classes of subroutines are named in the 
same manner as a FORTRAN variable (see the sec- 
tion, "Variables"). 

1. A subroutine name consists of 1 - 6 alpha- 
meric characters, the first of which must be 
alphabetic. 

2. The type of the function, which determines 
the type of the result, may be defined as 
follows: 

a. The type of an arithmetic statement func- 
tion may be indicated by the name of the 
function or by placing the name in a type 
statement. 

b. The type of a FUNCTION subprogram 
may be indicated by the name of the 
function (if it is real or integer) or by 
writing the type (REAL, INTEGER, COM- 
PLEX, DOUBLE PRECISION, LOGICAL) 
preceding the word FUNCTION. In the 
latter case, the type, implied by name, 

is overridden. 

c. The type of a built-in function is indi- 
cated within the FORTRAN Processor 
and need not appear in a type statement 
(see columns 5 and 6 of Figure 6-1). 

3. The type of a SUBROUTINE subprogram is 
unimportant and need not be defined, since 
the type of results returned is dependent only 



on the type of the variable names in the dum- 
my argument list. 

DEFINING SUBROUTINES 

The method of defining each class of subroutines is 
discussed below. 

Arithmetic Statement Functions 

Arithmetic statement functions are defined by a 
single arithmetic statement and apply only to the 
source program containing the definition. 



General Form 



a = b 
where: 

1 . the letter a is a function name followed by parentheses 
enclosing its arguments, which must be distinct, non- 
subscripted variables, separated by commas. 

2. the letter b is an expression that does not involve sub- 
scripted variables. Any arithmetic statement function 
appearing in b must have been previously defined. 



Examples; 

FIRST (X) = A*X+B 
JOB (X,B) = C*X+B 
THIRD F(D) = FIRST (E)/D 
MAX (A, I) = A**I-B - C 
LOGFCT (A, C) = A**2.GE. C/D 

1. As many as desired of the variables appear- 
ing in b may be stated in a as the arguments 
of the function. Since the arguments are dum- 
my variables, their names, which indicate 
the type of the variable, may be the same as 
names appearing elsewhere in the program 

of the same type. 

2. Those variables included in b that are not 
stated as arguments are the parameters of 

. the function. They are ordinary variables. 

3. All arithmetic statement function definitions 
must precede the first executable statement 
of the source program. 

4. The type of each argument must be defined 
preceding its use in the arithmetic state- 
ment function definition. 

Built-in Functions 

Built-in functions are pre-defined, open subroutines 
that exist within the FORTRAN Processor. A list 
of all the available built-in functions is given in Fig- 
ure 6-1. 



22 







Number of 




Type of 


Function 


Definition 


Arguments 


Name 


Argument 


Function 


Absolute value 


lArgI 


1 


ABS 


Real 


Real 








lABS 


Integer 


Integer 


Truncation 


Sign of Arg 


1 


AINT 


Real 


Real 




times largest 




INT 


Real 


Integer 




integer ^ |Arg | 










Remaindering 


Arg J (mod Arg2) 


2 


AMOD 


Real 


Real 


(see note below) 






MOD 


Integer 


Integer 


Choosing 


Max(Argj, 


>2 


AMAXO 


Integer 


Real 


largest value 


Arg2, . . • ) 




AMAXl 


Real 


Real 








MAXO 


Integer 


Integer 








MAXl 


Real 


Integer 


Choosing 


Min(Argi, 


>2 


AMINO 


Integer 


Real 


smallest value 


Arg2, • • • ) 




AMINl 


Real 


Real 








MINO 


Integer 


Integer 








MINI 


Real 


Integer 


Float 


Conversion from 
integer to real 


1 


FLOAT 


Integer 


Real 


Fix 


Conversion from 
real to integer 


1 


IFIX 


Real 


Integer 


Transfer 


Sign of Arg 2 


2 


SIGN 


Real 


Real 


of sign 


times 1 Argj | 




ISIGN 


Integer 


Integer 


Positive 


Arg J - Min 


2 


DIM 


Real 


Real 


difference 


(Arg^, Arg2) 




IDIM 


Integer 


Integer 


Obtain most 




1 


SNGL 


Double 


Real 


significant part 












of double - 












precision 












argument 












Obtain real 




1 


REAL 


Complex 


Real 


part of complex 












argument 












Obtain 




1 


AIMAG 


Complex 


Real 


imaginary part 












of complex 












argument 












Absolute value 


|Arg| 


1 


DABS 


Double 


Double 


Truncation 


Sign of Arg 
times largest 
integer < [Arg 


1 


IDINT 


Double 


Integer 


Choosing 


Max (Argp 


>2 


DMAXl 


Double 


Double 


largest 


Argg, . . . ) 










value 












Choosing 


Min (Arg J, 


>2 


DMINl 


Double 


Double 


smallest 


Arg2>...) 










value 












Transfer 


Sign of Arg 2 


2 


DSIGN 


Double 


Double 


of sign 


times 1 Arg J | 










Express single- 


D=(Arg, 0) 


1 


DBLE 


Real 


Double 


precision 












argument 












in double- 












precision form 












Express two real 


C=Arg]^+iArg2 


2 


CMPLX 


Real 


Complex 


arguments in 












complex form 












Obtain conjugate 


For Arg=X+iY, 


1 


CONJG 


Complex 


Complex 


of a complex 


C=X-iY 










argument 













NOTE: The function MOD (Arg J, Arg2) is defined as Arg^^ - lArgj/Arg J Arg2, where [[x] 
is the integral part of x 

Figure 6-1- Built-in Functions 



Subroutines, Functions, and Subprogram Statements 23 



FUNCTION Subprogram 

FUNCTION subprograms are defined by a special 
FORTRAN Source language program. 



General Form 



FUNCTION name (ai, a2, . . . , an) 

REAL FUNCTION name (aj, a2, . . , an) 

INTEGER FUNCTION name (a^, a2, . , , , an) 

DOUBLE PRECISION FUNCTION name (a^, a2, . . . , an) 

COMPLEX RTNCTION name (aj, ag, . . . , an) 

LOGICAL FUNCTION name (a^, a2, . . . , an) 

where: 

1. name is the symbolic name of a single-valued 
function, 

2. the arguments aj, a2, . . • i an, of which there must be at 
least one, are non- subscripted variable names or the 
dummy name of a SUBROUTINE or FUNCTION subprogram, 
and 

3. the type of the function may be explicitly stated preceding 
the word FUNCTION, 



Examples: 

FUNCTION ARCSIN(RADIAN) 

REAL FUNCTION ROOT (A,B,C) 

INTEGER FUNCTION CONST (ING, SG) 

DOUBLE PRECISION FUNCTION DBLPRE(R, S, T) 

COMPLEX FUNCTION CCOT(ABI) 

LOGICAL FUNCTION IFTRU (D, E, F) 



4. When a dummy argument is an array name, 

a DIMENSION or COMMON (with dimensions) 
statement must appear in the FUNCTION sub- 
program; also, the corresponding actual ar- 
gument must be a dimensioned array name. 

5. None of the dummy arguments may appear in 
an EQUIVALENCE statement in the FUNC- 
TION subprogram. 

6. The FUNCTION subprogram must be logic- 
ally terminated by a RETURN statement (see 
the section, "The RETURN Statement"). 

7. The FUNCTION subprogram may contain any 
FORTRAN statements except SUBROUTINE 
or another FUNCTION statement. 

8. The arguments of a FUNCTION subprogram 
may be any of the following: 

a. Any type of constant. 

b. Any type of subscripted or non-subscripted 
variable. 

c. An arithmetic or a logical expression. 

d. The name of a FUNCTION or SUBROU- 
TINE subprogram. 

9. A FUNCTION subprogram is referenced by 
using its name as an operand in an arithmetic 
expression. 

Those FUNCTION subprograms that are 
available with FORTRAN are given in Fig- 
ure 6-2. 

SUBROUTINE Subprogram 



1. The FUNCTION statement must be the first 
statement of a FUNCTION subprogram. 

2. The name of the function must appear at 
least once as a variable on the left side of an 
arithmetic statement or in an input statement. 

For example: 

FUNCTION CALC (A, B) 



CALC=Z+B 



RETURN 

By this means the output value of the function 
is returned to the calling program. 
The arguments may be considered dummy 
variable names that are replaced at the time 
of execution by the actual arguments supplied 
in the function reference in the calling pro- 
gram.. The actual arguments must corre- 
spond in number, order, and type with the 
dummy arguments. 



SUBROUTINE subprograms are defined by a special 
FORTRAN source language program. 



General Form 



SUBROUTINE name (ai, a2, . . . , an) 
where: 

1 . name is the symbolic name of a subprogram; and 

2. each argument, if any, is a non-subscripted variable 
name or the dummy name of a SUBROUTINE or FUNCTION 
subprogram.. 



Examples: 

SUBROUTINE MATMPY (A, N, M, B, L, J) 
SUBROUTINE QDRTIC (B,A,C, ROOTl,ROOT2) 

1. The SUBROUTINE statement must be the 
first statement of a SUBROUTINE subprogram. 

2. The SUBROUTINE subprogram may use one 
or more of its arguments to return output. 
The arguments so used must appear on the 
left side of an arithmetic statement or in an 
input list within the subprogram. 

3. The arguments may be considered dummy 
variable names that are replaced at the time 
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of execution b}' the actuai arguments suDplied 
in the CALL statement, which refers to the 
SUBROUTINE subprograin. The actual ar- 
guments must correspond in number, order, 
and type with the dummy arguments. 

Wlien a dummy argument is an array vnunn, 
a DLMIONSION or COMMON (with dimensions) 
statement must appear in the SUBROUTTNrC 
subprogram; also, the corresponding nctual 
argument in the CALL statement must be a 
dimensioned, array name. 

None of the dummy arguments may appear in 
an EQUIVALENCE statement in the SUBROU- 
TINE subprogram. 



6. The SUBROUTINE subprogi^am must Idc logi- 
cally terminated by a RETURN statcrftent. 

7. The SUBROUTINE subprogram may contain 
any FORTRAN,statomcnts except FUNCTION 
or another SUBROUTINE statement. 

Subprogram Names as Arguments 

FUNCTION and SUBROUTINE subprogram names 
may be the actual arguments of subprograms. In 
order to distinguish these subprogram names from 
ordinary variables when they appear in an argument 
list, they must appear in an EXTERNAL statement. 

EXTERNAL SIN 
CALLSUBR (A.SIN, B) 







Number of 




Type 


of 


Function 


Definition 


Arguments 


Name 


Argument 


Function 


Exponential 


e^rg 




EXP 


■Real 


Real 


Natural 


logg(Arg) 




ALOG 


Real 


Real 


logarithm 












Common 


logjQ(Arg) 




ALOGIO 


Real 


Real 


logaritlim 












Arctangent 


arctan (Arg) 




ATAN 


Real 


Real 




arctan (Arg^/ 




ATAN2 


Real 


Real 




Arg 2) 










Trigonometric 


sin (Arg) 




SIN 


Real 


Real 


sine 












Trigonometric 


cos (Arg) 




COS 


Real 


Real 


cosine 












Hyperbolic 


tanh (Arg) 




TAN'H 


Real 


Real 


tangent 


1/2 
(Arg) 










Square root 




SORT 


Real 


Real 


R c ill ainde ring 


Arg| {modArg2) 




DMOD 


Double 


Double 


Exponential 


eArg 




DEXP 


Double 


Double 


Natural 


logg (Arg) 




DLOG 


Double 


Double 


logaritlim 












Common 


logio(Arg) 




DLOGIO 


Double 


Double 


logarithm 












Arctangent 


arctan (Arg) 




DATAN 


Double 


Double 




arctan (Arg^/ 




DATAN2 


Double 


Double 




Arg 2) 










Trigonometric 


sin (Arg) 




DSIN 


Double 


Double 


sine 












Trigonometric 


cos (Arg) 




DCOS 


Double 


Double 


cosine 
Square root 


For Arg=X+iY 




DSQRT 


Double 


Double 


Absolute value 


(xKy^)'^^ 




CABS 


Complex 


Complex 


Exponential 


eArg 




CEXP 


Complex 


Complex 


Natural '''- 


logg(Arg) 




CLOG 


Complex 


Complex 


logarithm 












Trigonometric 


sin (Arg) 




CSIN 


Complex 


Complex 


sine 












Trigonometric 


cos(Arg) 




CCOS 


Complex 


Complex 


cosine 












Square root 


(Arg)^/2 




CSQRT 


Complex 


Complex 



Figure 6-2. Mathematical Subroutines 
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T U E C A L L S IW T E IM EN T 

The CALL statemont is used to refer to a SUBROU 
TINE subprog-ram. 



General Forni 



CALL SUBR (ap 32, . . , a,J 
where: 

1. SUBR is the name of a SUBROUTLNE subprogram, and 

2. aj,a2, ..,ajj are the n arguments. 



Examples: 

CALLMATMPY (X, 5, 10, Y, 7, 2) 

CALL QDRTIC (9. 732, Q/4. 536, R-S**2. 0, XI 

X2) 



Mathematical Subroutines 

! i'ORTRAN provides various commonly used mathe- 
j niatic;!! subroutines; these are defined as FUNCTION 
I subprog)\ams. The names of all of these subprograms 
j :ire automaticaily typed by the FORTRAN IV Compiler; 
! therefore, they need not appear in Type statements. 
I v^i-i-''-<^'^^''^s used as arguments of mathematical subrou- 
tines must be typed, either explicity or implicitly, in 
uccurciance with the function in which they appear. 
The mathematical subroutines are listed in Figure 6-2. 

A h ichine Indicator Tests 

la the following list of machine indicator test sub- 
routines, assume that i is an integer expression 
nnd that j is an integer variable. These subroutines 
:vro. rcvforenccd by CALL statements. 



The CALL statement transfers control to the sub- 
program and presents it with the actual arguments. 
The arguments may bo any of the following: 

1. Any type of constant. 

2. Any type of subscripted or non-subscripted 
variable. 

3. An arithmetic or a logical expression. 

4. Alphameric characters. Such arguments 
must be preceded by nil, where n is the count 
of characters included in the argument, e.g. , 
9HEND POINT. Note that blank spaces and 
special characters are considered in the 
character count when used in alphameric 
fields .. 

5. The name of a FUNCTION or SUBROUTINE 
subprogram. 

The ^arguments presented by the CALL statement 
must agree in. number, order, type, and array size 
(except as explained under the DIMENSION state- 
ment) with the corresponding arguments in the SUB- 
ROUTINE statement of the called subprogram. 



S LITE (i): If i = 0, all sense lights will be turned 
off. If i =1,2,3, or 4, the corresponding sense 

light will be turned on. 

S LITlr jT (i, j): Sense light i will be tested and turn- 
ed off. The variable j will be set to 1 if i was on, 
or j will be set to 2 if i was off. 

SSVVTCH (l,j): Sense switch i is tested and j is set 
to 1 if i was down and j is set to 2 if i was up. 

QVERFL (j): j is set to 1 if a floating point overflow 
condition exists, or j is set to 2 if no overflow con- 
dition exists. The machine is left in a no overflow 

condition. 

DVCHK (j)j If the divide check indicator is on, j is 
sot to 1 and the divide check indicator is turned off; 
if the divide check indicator is off, j is set to 2. 

EXIT, DUMP, and PDUMP 



SUBPROGRAMS PROVIDED BY FORTRAN 

FORTRAN includes several commionly used subrou- 
tines that are available to the programmer. The 
mathematical subroutines that are provided arc de- 
fined as FUNCTION subprograms; the subroutines 
provided to test the status of the machine nidicators 
(the sense switches and the sense lights) are de- 
fined as SUBROUTINE subprograms. In addition, 
FORTRAN includes the SUBROUTINE subprograms 
EXIT, DUMP, and PDUMP.. "EXIT terminates job 
execution; DUMP dumps core storage and then ter- 
minates job execution; PDUMP dum.ps core storage 
and then continues execution. 



EXIT 

A CALL to the EXIT subprogram terminates the 
execution of any program by returning control to the 

.Monitor. 

DUMP 

A CALL to the DUMP subprogram by the statement 

CALL DUMP (Ai, Bi, Fi An, Bn, Fj,) 

causes the indicated limits of core storage to be 
dumped and execution to be terminated by returning 
control to the Monitor. 

1. .1 and B are variable data names that indi- 
cate the limits of core storage to be dumped; 
either A or B may represent upper or lower 
limits. 
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Y\ is an integer indicating the clump lommt 
desired: 

F - dump in octal 

1 dump as real 

2 dump as integer 

3 dump in octal with mnemonics 

If no arguments are given, all of core slcra 
is dumped in octal- 

If the last argument Fn is omitted, it is as- 
sumed to be equal to and the dunap will be 
octal. 



PDUMP 

A CALL to the PDUMP subprogram by the statement 

CALL PDUMP (Ai, Bi, Fi, . . . , A^, B,^ F,,) 
causes the indicated limits of core storage to be 
dumped and execution to be continued. The PDUMP 
arguments are the same as the DUMP arguments. 

THE BLOCK DATA SUBPROGRAM 

In order to enter data into a COMMON block during 
compilation, a special subprogram must be written. 
This special subprogram contains only the DATA, 
COMMON, DIMENSION, and Type statements asso- 
ciated with the data being defined. Data may be 
entered into labeled, but not unlabeled, COMMON 
by the BLOCK DATA subprogram. 
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General Form 



BLOCK DATA 



The BLOCK DATA subprogram may not con- 
tain any executable statem.ents. 
The first statement of this subprogram must 
be the BLOCK DATA statement. 
All elements of a COMMON block must be 
listed in the COMMON statement even though 
they do not all appear in the DATA statement; 
for example, the variable A in the COMMON 
statement in Figure 6-3 does not appear in 
the DATA statement. 
Data may be entered into more than one 
COMMON block in a single BLOCK DATA 
subprogram. 



BLOCK DATA 

COMMON /ELN/C,A,B/RMG/Z,Y 

DIMENSION B( 4), Z( 3) 

DOUBLE PRECISION Z 

COMPLEX C 

DATA (B(I), 1=1, 4)/l. 1,1. 2, 2*1. 3/, 0/(2.4,3.769)/, 

Z(1)/7.6498085D0/ 
END 



Figure 6-3 
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CHAPTER 7: THE SPECIFICATION STATEMENTS 



The specification statements perform the followino,- 
functions: they provide information about storage 
allocation and about the constants and variables nsvd 
in the program. 

The DIMENSION Statement 



General Form 



DIMENSION vi(ii), V2(i2), V3(i3), ... j 

where: | 

1. each V]-n is a subscripted variable, and ! 

2. each ij, is composed of 1, 2, or 3 unsigned integer | 
constants and/or integer variables (i may be a variable onlyj 
when the DIMENSION statement "appears in a subprogram). 



Examples: 

DIMENSION A(10), B(5,15), C(L,M) 
DIMENSION S(10), K(5,5,5), G(IOO) 
The DIMENSION statement provides the informa- 
tion necessary to allocate storage for arrays in the 
object program. The DIMENSION statement defines 
■ the maximum size of arrays, 

1. Each variable that appears in subscripted form 
in the source program must appear in a 
DIMENSION statement contained in the source 
program; however, if the dimension informa- 
tion for a variable is included in a COMjMON 
or a Type statement in the source program, it 
must not be included in a DIMENSION statement, 

2. The DIMENSION statement must precede the 
first appearance of each subscripted variable 
in an executable or DATA statement for which 
it specified the size, 

3. A single DIMENSION statement- may specify 
tlie dimensions of any number of arrays, 

4. Dimensions specified in a COMMON or a Type 
statement are subject to all the rules for the 
DIMENSION statement. ■ 

Adjustable Dimensions 

The name of an array and the constants that are its 
dimensions may be passed as arguments in a sub- 
program call. In this way, a subprogram may per- 
form calculations on arrays whose size is not deter- 
mined until the subprogram is called. Figure 7-1 
illustrates the use of adjustable dim.ensions. 

1. Variables may be used as dimensions of an 
array only in the DIMENSION statement of a 
FUNCTION or SUBROUTINE subprogram. Fcr 
any such array, the array name and all the 
variables used as dimensions must appear as 
arguments in the FUNCTION or SUBROUTINJ-; 
statement. 



! 

SUBROUTINE MAYMY (. . 


.,r,'l,m,...) 


DIMENSION. ..,R(L,M),. 




DO 100 1=1, L 





"igiire 7-1 

2. The adjustable dimensions may not be altered 
within the subprogram. 

3. The absolute dimensions must be specified 

in a DIMENSION statement of the calling pro- 
gram. 

4. The calling program passes the specific di- 
mensions to the subprogram. These specific 
dimensions are those that appear in the DI- 
MENSION statement of the calling program. 
Variable dimension size may be passed 
through more than one level of subprogram. 



The COMMON Statement 



Cpneral Form 


COMMON a, b, c, . . , /r/d, e, f, . . . /s/g, h, . . . 




where: 




1 . a, b, . . . are variables that may be dimensioned, 




and 




2. /r/ , /s/, . . , are variables that are block names. 





Examples: 

COMMON A , B , C/X/Q , R/YY/M , P , Q 
COMMON /Z/G, H, J//D, F 

Variables, including array names, appearing in a 
COMMON statement are assigned locations relative 
to the beginning of a particular common block. This 
COMMON area may be shared by a program and its, 
subprograms . 

1, If the variables appearing in a COMMON state- 
ment contain dimension information, they must not 
appear in a DIMENSION statement or in a' Type state- 
ment that contains dimension information. 

2, The locations in the COMMON area are as- 
signed in the sequence in which the variables appear 
in the COMMON statement, beginning with the first: 
COMMON statement of the program. 

3, Elements placed in COMMON may be placed 
in sct'iaratc blocks. These separate blocks may 
share space in core storage at object time. Blocks 
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arc given names and those with the same name oc- 
cupy the same space. 

4. COMMON Block Names. The symbolic name 
of a block, which is 1-6 alpliajiieric characters tlic 
fii'st of which is alphabetic, precedes the varlnhle 
names bclong'ing to the block. The block name is 
always embedded in slashes, e.g. , /BB/, It must 
not be the Svame as the name of any other subroutine 
which is part of the same job. There are two types 
of COMMON blocks: blank and labeled. 

a. Blank COMMON is indicated either by omit- 
ting the block name if it appears at the be- 
ginning of the COMMON statement, or by 
preceding the blank COMMON variable by 
two consecutive slashes. 

b. Labeled COMMON is indicated by preceding 
the labeled COMMON variables by the block 
name embedded in slashes. 

5. The field of entries pertaining to a block name 
ends with a new block nam^c or v/ith the end of the 
COMMON statement or with a blank COMMON des- 
ignation. 

6. Block name entries are cumulative throughout 
the program. For example, the COMMON state- 
ments 

COMMON A , B , C /R/D , E/S/F 
COMMON G,H/R/I/S/P 
have the s'ame effect as the statement 

COMMON A, B, C, G, H/R/D, E, I/S/F, P 

7. Blank COMMON may be any length. Lnbeied 
COMMON must conform to the foUov/ing sl/.e .;•;;- 
quiremenl: all COMMON blocks of a given nmne 
must have the same length in all the programis that 
arc executed together. 

8. Variables brought into a COMMON block 
through EQUIVALENCE statements may increase 
the size of the block (see the section, "The EquiYa- 
lence Statement"), 

9. Two variables in COMMON may not be made 
equivalent to each other, directly or indirectly. 

10. A double-word variable in COMMON must 
be placed such that its high-order part is an even 
number of words away from the first clement in 
COMMON. . 



The EQUIVALENCE Statement 



General Form 



EQUIVALENCE (a, b, c, . . , ), (d, e, f, ...),... 

where: 

a, b, c, d, e, f, . . . are variables tliat may be subscripted; 

these subscripts must be integer constants. The number 

of subscripts appended to a variable must be equal to the 

number of dimensions of the variable. 



Examples: 

DIMENSION B(5) , C(10,10), D(5,10,15) 
EQUIVALENCE (A,B(1), C(5,4)), (D(1,4,3),E) 

The EQUIVALENCP] statement controls the alloca- 
tion of data storage by causing two or more variables 
to share the same core storage location. 

1. An EQUIVALENCE statement may be placed 
anywhere in the source program. Each pair of paren- 
theses in the statement list encloses the names of 
two or more variables that are to be stored in the 
same location during execution of the object program; 
any number of equivalences (i.e. , sets of parentheses) 
may be given, 

2. In the preceding example, the EQUIVALENCE 
statement indicates that A, and the B and C arrays 
are to be assigned storage locations so that the ele- 
naents A, B(l), and C(5,4) are to occupy the same 
location. In addition, it also specifies that i)(l,4,3) 
and E are to share the same location. 

3. Quantities or arrays that are not mentioned in 
an EQUIVALENCE statement will be assigned unique 
locations. 

4. Locations can be shared only among variables, 
not among constants. 

5. The sharing of storage locations requires a 
knowledge of which FORTRAN statements will cause 
a new value to be stored in a location. There are 
four such statemients: 

a. Execution of an arithmetic statement stores 
a new value in the variable on the left side 
of the equal sign. 

b. Execution of an ASSIGN i TO n statement 
stores a new value in n. 

c. Execution of a DO statement or an implied 
DO in an input/output list sometimes stores 
a new indexing value . 

d. Execution of a READ statement stores new 
values in the variables mentioned in the 
input list. 

6. Variables brought into a COMMON block through 
EQUIVALENCE statements may increase the size of 
the iDlock indicated by the COMMON statements, as in 
the following example: 

COMMON /X/A,B,C 

DIMENSION D(3) 

EQUIVALENCE (B,D(1)) 
The layout of core storage indicated by this ex- 
ample (extending from the lowest location of the block 
to the highest location of the block) is: 

A ' . 

B, D(l) 

C, D(2) 
D(3) 

7. Since arrays must be stored in consecutive for- 
ward locations , a variable may not be made equivalent 
to an element of an array in such a way as to cause 
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the array to exLcncl beyond the beginning of the 
COMMON block. 'For example, the following coclin.L: 
is invalid: 

COMMON /X/A,B,C 

DIMENSION D(3) 

EQUI\^7\LENCE (B,D(3)) 
l.)ecausc it would force D(l) to precede A, as foUuv-,-'-.: 

A, D(2) 

B, D(3) 
C 

8. The rule for making double-word variables 
equivalent to single-word variables is: 

In COMMON, the effect of Lhe EQUIVALENCE 
statements must be such that the high-order word of 
any doulDle-word variable is an even number of loca- 
tions away from the start of the COMMON block. 

In non-COMMON, the effect of the EQUlYALENCi'. 
statements must be such that the high-order word oi' 
any double-word variable is an even number of words 
away from the start of any other double-word varial)ic 
linked to it through EQUIVALENCE statements. 

9. Two variables in one COMMON block or in t^xu 
different COMMON blocks may not be made equiv- 
alent. 

The Type Statements 

The type of a variable or function may be specified 
by means of one of the six Type statements: 



General Form 



INTEGER a{ii), b(i2), 0(13),... 

REALa{ii), b(i2), c(i3), ... 

DOUBLE PRECISION a(ii), b{i2), c(i3), . . . 

COMPLEX a(ii), b{i 2), 0(13),.,. 

LOGICAL a(ii), b(i2), c{i^),... 

EXTERNAL a, b, c, . . . 

where: 

1. a, b, c, • . . are variable or function names appearing within 
the program. 

2. each i^ is comp>osed of 1, 2, or 3 integer constants and/or 
integer variables. Subscripts may only be appended to 
variable names appearing within the program, not function 
names. 



Examples: 

INTEGER BIXF, X, QF, LSL 

REALIMIN, LOG, GRN, KLW 

DOUBLE PRECISION Q, J, DSIN 

EXTERNAL SIN, MATMPY, INVTRY 

INTEGER A(10, 10), B 

COMPLEX C(4,5,3), D ■ 
The variable or function names following the type 
(INTEGER, REAL, etc.) in the Type statement. are 
defined to be of that type , and remain that type 



ihroughout the program; the type may not be 

changed. '-> '-.z^.'-r. 

N(;lc ihat LSL and GRN need not appear in their 
respective Type statements, , since their type is im- 
jjlicd Ijv their first characters. Note also that DSIN 
need not appear in its statement if it is used as a 
function in the program, since mathematical sub- 
rouimes are automatically typed by the FORTRAN IV 
Compiler, 

1, The appearance of a name in any Type state- 
jnent, except EXTERNAL, overrides the implicit 
ty pe a s s i gn m e nt . 

2, Variables that appear in EXTERNAL state- 
ments are subprogram names that are arguments of 
other suljprograms. 

3, A name may appear in two Type statements 
only if one of them is EXTERNAL. 

4, The Type statement must precede the first 
use u'i a name in any executable statement or DATA 
statement in the program. 

5, A name declared to be of a given type may 
iissurne onh' the values of a constant of the same type. 

G. The EXTERNAL statement may not be used 
to dimension variables. 

7, i\ny variable that is dimensioned by a Type 
statement may not be dimensioned elsewhere, i.e. , 
it may not appear in a DIMENSION statement or in 
a COMMON statement that contains dimension infor- 
mation. 

The DATA Statement 

Data may be compiled into the object program by 
means of the DATA statement. 



General Form 



DATA list/di,d2, ...,dn/,list/di,d2,k*d3, ...,dni/,; 
where: 

1. list contains the names of the variables being 
defined, 

2. d is tlie information literal, and 

3. k is an integer constant. 



]\xamples: 

DATA R, Q/14. 2 , 3HEND/, Z/0777777700001/ 
DATA (B(I),C(I), 1=1, 40, 2)/2. 0,3. 0,38*100.0/ 
LOGICAL LA, LB, LC, LD 
DATA LA, LB, LC, LD/F, . TRUE. ,. FALSE. , T/ 

1. List. Subscripted variables may appear in 
the list. Where a subscript symbol is used, it must 
be under control of DO-implying parentheses and 
associated parameters. Subscripts not so controlled 
must be integer constants. The DO-defining param- 
eters must be integer constants. 

2. jc. k may appear before a d-field to indicate 
that the field is to be repeated k times, k must be 



30 



Form C 28 -6 274-1 
Page Revised 11/29/63 
by TNL N28-0069 



followed ])y an * to separate it from the field to be 
repeated, 

3. d_. The d-literals may take any of the four 
following forms: 

a. Integer, real, double-precision, and eom- 
plex constants. 

b.' Alphameric characters. The alphameric 
field is written as nPI followed by n alpha- 
meric characters. Each group of six alpha- 
meric characters forms a word. If n is 
■ ,not a multiple of six, the remaining char- 
acters are left justified in the word, and 
the word is filled out with BCD blanks. 
Blanks are significant in alphameric fields. 

c. Octal digits. The octal field is written as 
O, followed by 1-12 signed or unsigned 
octal digits. 

d. Logical constants. The logical field may be 
written as either .TRUE. , . FALSE. , T, 
or F. 

4, There must be a one-to-one correspondence 
between the list items and the data literals. Each 
data literal (integer constant, real constant, alpha- 
meric constant, complex constant, logical constant, 
double-precision constant, or octal constant) corre- 



sponds to one element (undimensioned variable or 
subscripted array reference). Caution: If it is de- 
sired to define IG alphameric characters, say 
IGHDATAbTObBEbREA'Db starting at 0(1), then G 
must be dimensioned to cover at least three loca- 
tions and the entire literal corresponds to G(l). 

5. The BLOCK DATA subprogram (see page 27), 
which includes a DATA statement, compiles data 
into the common area of the program. 

fi. When DATA defined variables are redefined 
during execution, these variables will assume their 
new values regardless of the DATA statement. 

7. Where data is to be compiled into an entire 
array, the name of the array (with indexing informa- 
tion omitted) can be placed in the list. The number 
of data literals must be equal to the size of the 
array. 

For example, the statements 

DIMENSION B(25) 

DATA A,B,C/24*4.0,3.0,2.0,1.0/ 
define the values of A, B(l), .'. . , B(23) to be 4. 0, and 
the values of B(24), B(25), and C to be 3.0, 2.0, and 
1,0, respectively, 

8. The DATA statement may not be used to enter 
data into unlabeled COMMON. 
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APPENDIX A: SOURCE PROGRAM STATErJEN'J'S 
AND SECJUENCING 

The order in which the source prograin state me ills 
of a FORTRAN prograra are executed follows these 
rules: 

1. Control orjo-inates at the first executable 
statement. The specification statements 
and the FORJ.TAT, FUNCTION., SUBROU- 
TINE, and BLOCK DATA statement, are 
non-executable; m questions of scquencinti;, 
Lhcy can be ignored. 

2. If control was with statement S, then concroi 
will pass to the statement indicated by the 
normal sequencing properties of S. If, how- 
ever, S is ihe last statement in the range of 
one or more DOs that are not yet satisfied, 
then the normal sequencing of S is ignored 
and DO-sequencing occurs. 

The normal sequencing properties of each FOR- 
TRAN statement follows. 



Statement 

a = b 

ASSIGN t TO n 
BACKSPACE i 
BLOCK DATA 
CALL 

COMMON 
COMPLEX 
CONTINUE 



Normal Sequencing 

Next executable staiement 

Next executable statement 

Next executable statement 

Non-executable 

First statement of called 

subprogram 

Non- executable 

Non-executable 

Next executable statement 



Statement 



DIMENSION 
DO 

DOUBLE PRECISION 

END 

END FILE 

EQUIVALENCE 

EXTERNAL 

)-'OK,viAT 

FUNCTION 

GO TO n 

GO TO n, (n^jHo, . .. ,nm) 

GO TO (ni,n2,.t.,nm),i 

ir (a) ni,n2,n3 



IF (t) s 



INTEGER 
LOGICAL 

PAUSE 

i'iaNT 

PUNCH 
RFAD 
REAL 
RETURN 



REWIND 
STOP 

SUBROUTINE 
WRITE 



Normal Sequencing 

Non-executable 

Non-executable 

DO-sequencing, then the next 

executable statement 

Non-executable 

Terminates program 

Next executable statement 

Non-executable 

Non-executable 

Non-executable 

Non-executable 

Statement n 

Statement last assigned to n 

Statement nj 

Statement nj^,n2,n3 if a < 0, 

a=0, or if a > 0, respectively 

Statement s or next executable 

statement if t is true or false, 

respectively 

Non-executable 

Non-executable 

Next executable statement 

Next executable statement 

Next executable statement 

Next executable statement 

Non-executable 

The first statement, or part of a 

statement, following the reference 

to this program 

Next executable statement 

Terminates the program 

Non- executable 

Next executable statement 
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APPENDIX B: TABLE OF SOURCE PROGRAM 
CHARACTERS 



^1 
52 o 



1 01 01 A 



u 




aj 


0) 




a 


-t-i 




03 


^ 




Eh 




TD 

^ 


P 


4:3 


a 











pq 






12 
1 






?H 0) ^H 0) 

oo^Suo^Soo 
pqS^OOpq^oOpqS 



11 
61 21 J 1 



41 41 / 







21 61 



2 2 02 02 B 



3 3 03 03 C 



12 
2 

12 
3 



11 

62 22 K 2 42 42 

11 



63 23 







22 62 



43 43 T 3 23 63 



4 4 04 04 D 



12 

4 



11 







64 24 M 4 44 44 U 4 24 64 



5 5 05 05 E 



12 11 

5 65 25 N 5 




45 45 V 5 25 65 



12 
6 6 06 06 F 6 66 26 O 



11 

6 46 46 W 6 26 66 



7 7 07 


07 


G 


12 

7 


67 


27 


P 


11 

7 


47 


47 


X 




7 


27 


67 


8 8 10 


10 


H 


12 

8 


70 


30 


Q 


11 

8 


50 


50 


Y 




8 


30 


70 


9 9 11 


11 


I 


12 
9 


71 


31 


R 


11 
9 


51 


51 


Z 



9 


31 


71 


blank blank 20 


60 


+ 


12 


60 


20 


- 


11 


40 


40 








12 


00 


= 8-3 13 


13 


, 


12 
8-3 


73 


33 


$ 


11 
8-3 


53 


53 


f 



8-3 


33 


73 


' 8-4 14 


14 


) 


12 

8-4 


74 


34 


* 


11 

8-4 


54 


54 


i 




8-4 


34 


74 


The characters $ and 
text in an H field. 


' (8-4) can be used 


in FORTRAN only as 


alphameric 





APPENDIX C: DIFFERENCES BETWEEN FORTRAN 
II and FORTRAN IV 

1. All langu^e items distinguished by a column 
1 modal punch, except B, in FORTRAN II have been 
incorporated into FORTRAN IV by the Type state- 
ments as follows: 



FORTRAN II Modal Punch 



FORTRAN IV Type Statement 
(See the section, "The 
Type Statements") 



Double-precision arithmetic - D DOUBLE PRECISION 



Complex arithmetic - I 



COMPLEX 



F-card - F 



EXTERNAL 
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M \/£ 



T 



The DATA atatement may be used to enter octal 
constants intc/a FORTRAN IV program; the pro- 
grammer fli^ptivrite subroutines, using these con- 
stants, to handle Boolean arithmetic. 

2. Function Naming. 

a. Where the initial character of a function 
name is used to denote the type as float- 
ing point (real) or fixed point (integer) in 
FORTRAN II, incompatibilities may arise. 
In FORTRAN IV, this difficulty is han- 
dled by the TYPE statements REAL and 
INTEGER, which define a variable name 
or function name as floating point or fixed 
point, respectively (see the section, "The 
Type Statements"). 

b. The number of characters in an open, a 
closed, or an arithmetic statement function 
name in FORTRAN II is 4 through 7 , ending 
in F; whereas, in FORTRAN IV, the num- 
ber of characters is 1 through 6 and the 
final F has no meaning. In both cases, the 
first character of the function name must 
be alphabetic (see chapter 6). 

c. Built-in and arithmetic statement functions 
are not identified by a terminal F in FOR- 
TRAN IV; they are named in FORTRAN 
rv as described in b. above. The FOR- 
TRAN II library function is a FORTRAN 
rv FUNCTION subprogram which is in- 
ternal to the Processor. 

3. COMMON and EQUIVALENCE. 

a. In FORTRAN IV, EQUIVALENCE does 
not effect the ordering within COMMON, 
and it does not create a gap in COMMON 
storage; the only effect it can have on a 
COMMON block is to make its size greater 
than that indicated by the COMMON 
statements of the program (see the sec- 
tion, "The COMMON Statement"). 

b. The FORTRAN IV COMMON statement 
may contain dimension information. 

4. In FORTRAN IV, if an explicit type is given to 
a variable name which is used throughout the pro- 
gram as an ordinary variable and also as a dummy 
argument of an arithmetic statement function, the ex- 
plicit type applies in both contexts. 

5. Implicit multiplication, which occurs in FOR- 
TRAN II as a by-product of the arithmetic translator 
techniques, is not permitted in FORTRAN IV. Thus, 
the following combinations are not permitted in FOR- 
TRAN IV: 

K( ) 

( )V 

( )K 
where V is a variable, K is a constant, and ( ) is 
any arithmetic expression within parentheses. 



6. The FORTRAN II statements in column 1 are 
changed to the FORTRAN IV statements in column 2. 



FORTRAN II Statement 



IF ACCUMULATOR OVERFLOW n, n2 



IF QUOTIENT OVERFLOW in,n2 



IF DIVIDE CHECK n, n2 



IF (SENSE SWITCH i) ni, x^2 



SENSE LIGHT i 



IF (SENSE LIGHT i) ni,n2 



READ TAPE i, list 



READ INPUT TAPE i,n,list 



WRITE TAPE i, list 



WRITE OUTPUT TAPE i, n, list 



FORTRAN IV Statement 



CALLOVERFL(j) 



CALL OVERFL (j) 



CALL DVCHK (j) 



CALLSSWTCH (i, j) 



CALL SLITE (i) 



CALLSLITET (i,j) 



READ (i) list 



READ (i,n) list 



WRITE (r) list 



WRITE (i.n) list 



Binary record 
BCD record 
Binary record 
BCD record 



The FREQUENCY, READ DRUM, WRITE DRUM statements 
of FORTRAN II are not part of the FORTRAN IV language. 

7. Additional FORTRAN IV statements. 

a. DATA (see the section, "The DATA 
Statement"). 

b. BLOCK DATA (see the section, "The 
Block Data Subprogram"). 

c. LOGICAL, an additional Type statement 
which defines variables to be used in logi- 
cal computation (see the section "The 
Type Statement"). 

8. Differences in Output Produced by FORTRAN 
IV and FORTRAN II. 

Programmers will find that the output produced by 
a source program in FORTRAN IV may not be the 
same as that provided by the identical program in 
FORTRAN II. Where differences do occur, they are 
attributed exclusively to the following differences 
which exist between FORTRAN IV and FORTRAN II. 

a. The logarithm subroutine of FORTRAN IV 
employs a new algorithm which yields more 
accurate results for most arguments than 
does the logarithm subroutine of FORTRAN 
II. 

b. Floating point constants which are written 
into the source program are converted by 
FORTRAN rv by a somewhat different 
algorithm than that used by FORTRAN IL 
The result is that more significance tends 
to be preserved and a more accurate con- 
version is achieved by FORTRAN IV as 
compared to its predecessor. 

c. The mathematical subroutines in FORTRAN 
IV are assembled by MAP, and those in 
FORTRAN II are assembled by FAP. The 
conversion routines in MAP provide more 
precise conversions for constants than do 
those in FAP. As a consequence, FOR- 
TRAN rv tends to produce more precise 
results than FORTRAN II for those sub- 
routines which use the same algorithm 
(and its associated constants) . The SIN/ 
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COS subroutine is a very good example of 
this effect, 
d. The order in which a sequence of multi- 
plications (or of multiplications and 
divisions) is executed by the object pro- 
gram in FORTRAN IV may be different 
from that in FORTRAN II. If such a dif- 
ference in ordering should occur, neither 
method may be considered superior to the 
other from the standpoint of computational 
accuracy. 



APPENDIX D: 
GRAM SIZE 



LIMITATIONS ON SOURCE PRO- 



In translating a source program into an object pro- 
gram, FORTRAN forms and utilizes various tables 
containing information about the source program. 
These tables are divided into two categories. 

First Category Tables 

The first category of tables consists of those tables 
for which space is allocated only when a table is 
actually required. 

Storage of First Category Tables 

At such time as a first category table is required, a 
buffer is assigned by the Table Routine in the TR 
buffer area, which begins immediately after the 
FORTRAN compiler in core storage and extends 
toward the top of storage. The top of the TR buffer 
area is determined by the NAME table. 

The NAME table is the only first category table 
that is not stored in the TR buffer area. Rather, 
the NAME table storage begins near the end of core 
storage and extends downward toward the TR buffer 
area. 

The following map of core storage indicates the 
location and direction of the TR buffer area and the 
NAME table. 



NAME 



TR Buffer Area 



FORTRAN Compiler 



77777 

Beginning of NAME Table 



Beginning of TR Buffer Area 



Essentially, the TR buffers and the NAME table 
share the same space but start filling the space 
from opposite ends. When they meet an overflow 
occurs. 

Overflow of TR Buffer Area and NAME Table Storage 

It is unlikely that the TR buffers and the NAME table 
would overflow each other; however, if they did, a 
failure to compile would result. 

If such an overflow occurs, it would be possible 
to recompile the program only after reducing the 
size of the TR buffers or by modifying the program 
so that excessive table entries will be eliminated. 

First Category Tables 

The following tables belong to the first category. 
With the following description of what causes entries 
to be made in the tables , the programmer should be 
able to reduce the size of tables if an overflow occurs. 

The following comments and definitions apply to 
the tables: 

1. The phrase "literal appearance" means that 
each item must be counted each time it 
appears . 

2. The BETA, LAMBDA, ASFD, and RANG E2 
tables are used for one statement only and 
are released immediately upon completion of 
processing of that statement. 

3. In computing the table size of BETA, LAMBDA, 
and RANGE 2 , all subscripts are ignored. 

4. A subexpression is any term or expression 
within parentheses; a function argument is a 
subexpression. 

5. Each of the FORTRAN operators is given a 
value when computing the size of the BETA 
and LAMBDA tables . The operator .NOT. 
is ignored. 



Operator 

.OR. 

.AND. 

Relational 

Operators 

+,- 
*,/ 



Value 

1 

2 



Function 

Reference, 7 

The symbol © for a function reference is used as 
follows: MAX (X, Y) in a FORTRAN statement is 
interpreted by the compiler as if it had been written 
MAX @ X © Y; of course, the programmer does 
not use the function symbol @ . 

6. [x] is the greatest integer less than or equal 
to X. 
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ASFD: A one word entry is made for each argument DELTA Table; An entry is made for each assigned 

and computed GO TO. The number of words in 



of an arithmetic statement function definition. 

BETA Table : The size of the BETA table is com- 
puted by the following procedure: 

1. Ignoring all subexpressions: 

a. Scan the right side of the arithmetic 
statement for the first operator; add the 
value of this operator to a counter. 

b. Scan for the next operator; if the value 
of this operator is less than or equal to 
the value of the immediately preceding 
operator, continue scanning; if the value 
of this operator is greater than the im- 
mediately preceding operator, add the 
difference to the counter; continue scan- 
ning and totaling until the entire expres- 
sion has been processed. 

2. For each subexpression not an argument, 
add to the counter the difference between 7 and the 
value of the operator which immediately precedes 
the left parenthesis of the subexpression. 

3. For each subexpression, previously ignored, 
perform the processing described under item 1. 
above, keeping a continuous total in the counter. 

4. The value in the counter is the table size. 



BR ADD Table : The number of entries is equal to 
the number of entries in the BRANCH table. 

BRANCH Table: An entry is made for each state- 
ment that may result in a transfer of control. The 
number of words in each entry is found in the fol- 
lowing table: 



Statement 



Number of Words in Entry 



Logical IF 
Arithmetic IF 
Unconditional GO TO 
Assigned GO TO 
Computed GO TO 



1 

3 

1 

number of statement numbers in list 

number of statement numbers in list 



COMMON Table; A one word entry is made for 
each literal appearance of variables in a COMMON 
statement. A one word entry is also made for each 
specification of blank or labeled COMMON. 

CONST Table; An entry is made for each literal 
appearance of a constant. The number of words in 
each entry is 2 +r-^]. where n is the number of 
non-blank characters used for the representation of 
the constant in the program. 

DDO Table; A two word entry is made for each im- 
plied DO in a DATA statement. 



each entry is 1+ (n/2) , where n is the number of 
statement numbers in the list. 

DIM Table; An entry is made for each assignment 
of a dimension to an array. The number of words 
in each entry is 1+ f^l , where n is the number of 
dimensions. 

DLIST Table; An entry is made for each DATA 
statement. The entry size is n + d + v + 2d where n 
is the number of lists in the DATA statement, v is 
the number of variables, and d is the number of 
implied DOs. 

DLIT Table; A one word entry is made for each 
list in a DATA statement. An entry is also made 
for each literal (constant) appearing in a DATA 
statement. The size of this latter entry is 2 + [o^^, 
where n is the number of non-blank characters in 
the literal. 

EQUIT Table; A one word entry is made for each 
literal appearance of variables in an EQUIVALENCE 
statement. 



FORMAT Table: 



An entry is made for each FOR- 
The number of words in each en- 



MAT statement, 
try is 2 + i^^l i where n is the number of charac- 
ters following the word FORMAT and blanks are re- 
tained in an alphameric field, but deleted elsewhere. 

FORT AG Table; A two word entry is made for eaich 
literal appearance of a subscripted variable in any 
statement except a DATA statement. 

FORVAL Table; A one word entry is made for each 
literal appearance of a non- subscripted integer vari- 
able on the left side of an arithmetic statement, in 
an I/O list, in COMMON statements, and in the 
argument list of a SUBROUTINE subprogram. 

INC DO Table; A one word entry is made for each 
left parenthesis that is not used for subscripts in a 
DATA statement. This table lasts for one statement 
only. 

LAMBDA Table: The size of the LAMBDA table is 
computed by the following procedure. 

1. Set a counter equal to the number of constants, 
variables , and functions in the arithmetic 
expression. 

2. Ignoring all subexpressions, scan the expres- 
sion for the first operator: if the next opera- 
tor-has the same value as the immediately 



36 



3. 



preceding operator, continue scanning; if the 
value of the next operator is less than the 
immediately preceding operator, add one to 
the counter and continue scanning; if the val- 
ue of the next operator is greater than the 
imm^ediately preceding operator, add the dif- 
ference in values to the counter and continue 
scanning and totaling until the entire expres- 
sion has been processed. 
For each subexpression, perform the pro- 
cessing described under item 2, keeping a 
cumulative value in the counter. 
The value in the counter is the table size. 



Second Category Tables 

The second category of tables consists of those 
tables for which space is specifically allocated in 
the Compiler program. The space allocated for 
each of these tables remains assigned even though 
the table may not be required. 

It is unlikely that a compilation failure will re- 
sult from an overflow of one of these tables. 

For a complete description of these tables and 
their limitations , see the FORTRAN section of the 
7090/7094 IB JOB Processor publication. 



NAME Table; A two word entry is made for each 
variable, function name, and subroutine name (other 
than the name of the program, being compiled) in the 
source program. 

RANGE2 Table; A one word entry is made for each 
function argument that is an expression. 

SUBARG Table; An entry is made for each FUNC- 
TION or SUBROUTINE subprogram. The number of 
words in each entry is 1 + [n/2J , where n is the 
number of subprogram arguments. 

TAU Table; A one word entry is m^ade for each ap- 
pearance of a unique subscript combination. Sub- 
script combinations that are actually the same as 
others are considered unique if there is a difference 
in the size of the dimensions of the subscripts. 

TDO Table; A three word entry is made for each DO 
or implied DO. 

TDOB Table; A one word entry is made for each 
DO or implied DO. 

TEIFNO Table; A one word entry is made for each 
executable statement which has a statement number. 



APPENDIX E:' OPTIMIZATION OF ARITHMETIC 
EXPRESSIONS 

In order to optimize the object program, a sequence 
of operations on the multiply -divide (*,/) level may 
be reordered. The reordering tends to an alterna- 
tion of the multiply and divide operations. It occurs 
where all elements of the expression are real type. 
This is done on the assumption that mathematically 
equivalent expressions are computationally equiva- 
lent. 

Where the multiply-divide expression involves 
naixed real and complex types, the operations on the 
real types occur first and are alternated. 

Where the order of operations is considered sig- 
nificant, the programmer may use nested parenthe- 
ses in the expression. 

There are two distinct ways a given floating-point 
constant may enter computation within a FORTRAN 
program. It may be 

1. written into the FORTRAN source program, or 

2. read in as data during execution. 

In the first version of the 7090/7094 FORTRAN 
Compiler (IBFTC) and its Library (IBLIB) , there 
may be a difference in the low-order binary bit of 
the same constant arising from these two sources. 
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the next sequential statement unless s is an 
arithmetic IF-or GO TO-type statement, in 
which case , control is transferred as indicated. 

2. If t is false, control is transferred to the next 
sequential statement. 

3. If t is true and s is a CALL statement, upon 
return from the subprogram control is trans- 
ferred to the next sequential statement. 

The DO Statement 



DO 



General Form 



DO n i = m< , m-, m3 

where: 

1 . n is a statement number, 

i is a non-subscripted integer variable, and 

3. mj , m2, m3 are each either an unsigned integer constant or a non- 
subscripted integer variable; if m3 is not stated, it istakentobe 1., 



Examples: 

DO 30 I = 1, M, 2 
DO 24 I = 1, 10 
The DO statement is a command to execute re- 
peatedly the statements that follow, up to and in- 
cluding the statement numbered n. The statements 
in the range of the DO are executed repeatedly with 
i equal to mi, then i equal to mi + ms, then i equal 
to mi + 2m3, etc. , until i is equal to the highest 
value in this sequence that does not exceed ui2. The 
statements in the range of the DO will be executed 
at least once. 

1. The range of a DO is that set of statements 
that will be executed repeatedly; i.e. , it is the se- 
quence of consecutive statements immediately follow 
ing the DO statement, up to and including the state- 
ment numbered n. After the last execution of the 
range, the DO is said to be satisfied. 

2. The index of a DO is the integer variable i. 
Throughout the range of the DO, the index is avail- 
able for computation, either as an ordinary integer 
variable or as the variable of a subscript. Upon 
exiting from a DO by satisfying the DO, the index i 
must be redefined before it is used in computation. 
Upon exiting from a DO by transferring out of the 
range of the DO, the index i is available for com- 
putation and is equal to the last value it attained. 

3. DOs within DOs . Among the statements in the 
range of a DO may be other DO statements; such a 
configuration is called a nest of DOs. If the range 
of a DO includes another DO, then all of the state- 
ments in the range of the latter must also be in the 
range of the former. 

4. Transfer of Control and DOs . Control may not 
be transferred into the range of a DO from outside 
its range. Thus, in the configuration below 1,2, and 
3 are permitted transfers, but 4,5, and 6 are not. 




a. 



b. 



Restrictions on Statements in the Range of a 
DO ~ 

Any statement that redefines the index or 
any of the indexing parameters (m's) is 
not permitted in the range of a DO. 
The range of a DO cannot end with an 
arithmetic IF- or GO TO-type statement, 
lie range of a DO may end with a logical 
IF, in which case, control is handled as 
follows: if the logical expression t is 
false, the DO is reiterated; if the logical 
expression t is true, statement s is exe- 
cuted and then the DO is reiterated. How- 
ever, if t is true and s is an arithmetic 
IF or transfer type statement, control is 
transferred as indicated. 
6. When a reference to a subprogram is executed 
in the range of a DO, care must be taken that the 
called subprogram does not alter the DO index or 
the indexing parameters . 

The CONTINUE Statement 



General Form 



CONTINUE 



CONTINUE is a dummy statement that gives rise 
to no instructions in the object program. It is most 
frequently used as the last statement in the range of 
a DO to provide a transfer address for IF and GO TO 
statements that are intended to begin another repe- 
tition of the DO range. 
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The PAUSE Statement 



General Form 



PAUSE or PAUSE n 

where n is an unsigned octal integer constant of 1-5 digits. 



Examples : 

PAUSE 

PAUSE 77777 
The machine will halt with the octal number n in 
the address field of the Storage Register. If n is not 
specified, it is understood to be zero. Depressing 
the START key causes the program to resume exe- 
cution of the object program with the next executable 
FORTRAN statement. 

The END Statement 



General Form 



END 



1. The END statement terminates compilation 
of a program. 



2. The END statement must be the physically 
last statement of the program. 

The STOP Statement 



General Form 



STOP 



The STOP statement terminates the execution of 
any program by returning control to the Monitor. 
Execution of a program may also be terminated by 
a CALL to the EXIT and DUMP subroutines (see 
the section, "EXIT, DUMP, and PDUMP"). 

The RETURN Statement 



General Form 



RETURN 



The normal exit from any subprogram is the 
RETURN statement, which returns control to the 
calling program. The RETURN statement is the 
logical end of the program; there may be any num- 
ber of RETURN statements in the program. 
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4. The field width w, for D-, E-, and F-conver- 
sion, must include a space for the decimal point and 

a space for the sign. Thus, for D- and E-conversion, 
w>d+7, and for F-conversion, w^d+S. 

5. The exponent, which may be used with D- and 
E-conversion, is the power of 10 to which the num- 
ber must be raised to obtain its true value. The ex- 
ponent is written with an E (for E-conversion) or D 
(for D-conversion) followed by a minus sign if the 
exponent is negative, or a plus sign or a blank if the 
exponent is positive, and then followed by two num- 
bers which are the exponent. For example, the num- 
ber .002 is equivalent to the number .2E-02. 

6. For input under D-conversion, up to 17 deci- 
mal digits are converted and the result is stored so 
that the most significant part and the least signifi- 
cant part are in adjacent core storage locations. 

For output under D-conversion, the two core stor- 
age words representing the double -precis ion quan- 
tity are considered one piece of data and converted 
as such. 

7. If a number converted by I-conversion re- 
quires more spaces than are allowed by the field 
width w, the excess on the high-order side is lost. 
If the number requires fewer than w spaces , the 
leftmost spaces are filled with blanks. If the number 
is negative , the space preceding the leftmost digit 
will contain a minus sign if sufficient spaces have 
been reserved. 

Complex Number Fields 



The basic difference between A- and H-conversion 
is that information handled by A-conversion is given 
a variable name or array name and hence can be re- 
ferred to by means of this name for processing and 
modification, whereas, information handled by H- 
conversion is not given a name and may not be re- 
ferred to or manipulated in storage in any way. 



A-Conversion 

The variable name to be converted by A-conversion 
must conform to the normal rules for naming FOR- 
TRAN variables; it may be any type. 

1. On input, nAw will be interpreted to mean 
that the next n successive fields of w charac- 
ters each are to be stored as BCD informa- 
tion. If w is greater than 6, only the 6 right- 
most characters will be significant. If w is 
less than 6, the characters will be left-ad- 
justed, and the word filled out with blanks. 

2. On output, nAw will be interpreted to mean 
that the next n successive fields of w char- 
acters each are to be the result of transmis- 
sion from storage without conversion. If w 
exceeds 6, only 6 characters of output will 
be transniitted, preceded by w-6 blanks. If 
w is less than 6, the w leftmost characters 
of the word will be transmitted. 



Since a complex quantity consists of two separate 
and independent real numbers, a complex number 
is transmitted by two successive real number speci- 
fications or by one real number specification that is 
repeated. 

See Figure 5-1 for an example of a FORMAT 
statement to transmit an array consisting of six 
complex numbers. 



FORMAT (2E10. 2, E8.3, 1PE9.4, E10.2, F8.4, 3(E10.2,F8.2)) 



Figure 5-1 

Alphameric Fields 

FORTRAN provides two ways by which alphameric 
information may be transmitted; both specifications 
result in storing the alphameric information inter- 
nally in BCD. 

1. The specification Aw causes w characters to 
be read into, or written from, a variable or 
array name. 

2. The specification nH introduces alphameric 
information into a FORMAT statement. 



H-Conversion 

The specification nH is followed in the FORMAT 
statement by n alphameric characters . For ex- 
ample, 

31 H THIS IS ALPHAMERIC INFORMATION 
Note that blanks are considered alphameric charac- 
ters and must be included as part of the count n. 
The effect of nH depends on whether it is used with 
input or output. 

1 . On input, n characters are extracted from 
the input record and replace the n characters 
included with the source program FORMAT 
specification. 

2. On output, the n characters following the 
specification, or the characters that replaced 
them, are written as part of the outfxit rec- 
ord. 

See Figure 5-2 for an example of A- and H-con- 
version in a FORMAT statement. 

The statement FORMAT (4HbXY-, F8.3, A8) might 
produce the following lines, where b indicates a 
blank character: 
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XY = b-93. 210bbbbbbbb 
XY = 9999.999bbOVFLOW 
XY = bb28.768bbbbbbbb 



Figure 5-2 



Figure 5-2 assumes that there are steps in the 
source program that read the data OVFLOW, store 
this data in the word to be printed in the format A8 
when overflow occurs, and store six blanks in the 
word jyhen overflow does not. occur. 

Blank Fields - X-Conversion 

The specification nX introduces n blank characters 
into an input/output record where 0:$n<132. 

1. On input, nX causes n characters in the input 
record to be skipped, regardless of what they 
actually are. 

2. On output, nX causes n blanks to be intro- 
duced into the output record. 

Repetition of F ield Format 

It may be desired to print or read n successive fields 
in the same format within one record. This may be 
specified by giving n, an unsigned integer, before 
D, E,F,I, O, or A. Thus, the field specification 
3E12. 4 is the same as writing El 2. 4, E12.4, E12.4. 



The scale factor is defined for output as follows: 
external quantity=intemal quantity 
X 10 sc^l6 factor 

For input, scale factors have effect only on F- 
conversion. For example, if input data is in the 
form xx.xxxx and it is desired to use it internally in 
the form .xxxxxx, then the FORMAT specification 
to effect this change is 2PF7.4. For output, scale 
factors may be used with D-, E-, and F-conversion. 

For example, the statement FORMAT (12, 3F11.3) 
might give the following printed line: 

2'|bbbb-93. 20gbbbbb-0. 008|DbbbbbO. 554 
But the statement FORMAT (12, 1P3F11. 3) used with 
the same data would give the following line: 
2|bbb-932. OO^fDbbbb-O. 076l)bbbbb5. 536 
Whereas, the statement FORMAT (12, -1P3F11.3) 
would give the following line: . 

21bbbbb-9. 321fcbbbb-0. OOlfbbbbbbO. 055 

A positive scale factor used for output with D-and 
E-conversion increases the number and decreases the 
exponent. Thus, FORMAT (12, 1P3E12.4) would 
produce with the same data; 

243-9. 3209Ebo4)-7.5804E-oA)b5.5536E-01 

The scale factor is assumed to oe zero if no other 
value has been given. However, once a value has 
been given, it will hold for all D-, E-, and F-con- 
versions following the scale factor within the same 
FORMAT statement. This applies to both single- 
record formats and multiple -record formats (see 
below). Once the scale factor has been given, a 
subsequent scale factor of zero in the same FOR- 
MAT statement must be specified by OP. Scale fac- 
tors have no effect on I- and O-conversion. 



Repetition of Groups 

A limited parenthetical expression is permitted in 
order to enable repetition of data fields according 
to certain format specifications within a longer FOR- 
MAT statement. Thus, FORMAT (2(F10.6,E10.2),I4) 
is equivalent to FORMAT (FIO. 6, ElO. 2, FIO. 6, 
EIO.2,14). Two levels of parentheses, in addition 
to the parentheses required by the FORMAT state- 
ment, are permitted. The second level of paren- 
theses facilitates the transmission of complex quan- 
tities. 

Scale Factors 

To permit more general use of D-, E-, and F-con- 
version, a scale factor followed by the letter P may 
precede the specification. The magnitude of the 
scale factor must be between -8 and +8, inclusive. 
The scale factor is defined for input as follows: 
10-scale factor x external quantity = internal 
quantity 



Multiple -Record Formats 

To deal with a block of more than one line of print, 
a FORMAT specification may have several different 
one-line formats separated by a slash, /, to indicate 
the beginning of a new blank line. Thus, FORMAT 
(3F9.2,2F10.4/8E14.5) would specify a multi-line 
block of print in which lines 1 , 3 , 5 , ... have format 
(3F9.2,2F10.4), and lines 2, 4, 6, .. . have format 
(8E14.5). 

If a multiple-line format is desired in which the 
first two lines are to be printed according to a spe- 
cial format and all remaining lines according to an- 
other format, the last line-specification should be 
enclosed in a second pair of parentheses; e.g. , 
FORMAT (I2,3E12.4/2F10.3,3F9.4/ (10F12.4)). 
If data items remain to be transmitted after the for- 
mat specification has been completely "used," the 
format repeats from the last previous parenthesis 
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5. 



of execution by the actual arguments supplied 
in the CALL statement, which refers to the 
SUBROUTINE subprogram. The actual ar- 
guments must correspond in number, order, 
and type with the dummy arguments. 

When a dummy argument is an array name , 
a DIMENSION or COMMON (with dimensions) 
statement must appear in the SUBROUTINE 
subprogram; also, the corresponding actual 
argument in the CALL statement must be a 
dimensioned, array name. 

None of the dummy argum.ents may appear in 
an EQUIVALENCE statement in the SUBROU- 
TINE subprogram. 



6. The SUBROUTINE subprogram must be logi- 
cally terminated by a RETURN statement. 

7. The SUBROUTINE subprogram may contain 
any FORTRAN statements except FUNCTION 
or another SUBROUTINE statement. 

Subprogram Names as Arguments 

FUNCTION and SUBROUTINE subprogram names 
may be the actual arguments of subprograms. In 
order to distinguish these subprogram names from 
ordinary variables when they appear in an argument 
list, they must appear in an EXTERNAL statement. 

EXTERNAL SIN 
CALL SUBR (A, SIN, B) 







Number of 




Type of 


Function 


Definition 


Arguments 


Name 


Argument 


Function 


Exponential 


eArg 




EXP 


Real 


Real 


Natural 


logg(Arg) 




ALOG 


Real 


Real 


logarithm 












Common 


logiO(Arg) 




ALOGIO 


Real 


Real 


logarithm 












Arctangent 


arctan (Arg) 




ATAN 


Real 


Real 




arctan (Arg^/ 




ATAN2 


Real 


Real 




Arg2) 










Trigonometric 


sin (Arg) 




SIN 


Real 


Real 


sine 












Trigonometric 


cos (Arg) 




COS 


Real 


Real 


cosine 












Hyperbolic 


tanh (Arg) 




TANH 


Real 


Real 


tangent 


(Arg) 










Square root 




SQRT 


Real 


Real 


Remaindering 


Argj^ (modArg2) 




DMOD 


Etouble 


Double 


Exponential 


gArg 




DEXP 


Double 


Double 


Natural 


logg (Arg) 




DLOG 


Double 


Double 


logarithm 












Common 


log 10 (Arg) 




DLOGIO 


Double 


Double 


logarithm 












Arctangent 


arctan (Arg) 




DATAN 


Double 


Double 




arctan (Argj/ 




DATAN2 


Double 


Double 




Arg 2) 










Trigonometric 


sin (Arg) 




DSIN 


Double 


Double 


sine 












Trigonometric 


cos (Arg) 




DCOS 


Double 


Double 


cosine 












Square root 


(Arg)^/2 

For Arg=X+iY 




DSQRT 


Double 


Double 


Absolute value 


(x2+y2)1/2 




CABS 


Complex 


Complex 


Exponential 


eArg 




CEXP 


Complex 


Complex 


Natural 


logg(Arg) 




CLOG 


Complex 


Complex 


logarithm 












Trigonometric 


sin (Arg) 




CSIN 


Complex 


Complex 


sine 












Trigonometric 


cos(Arg) 




CCOS 


Complex 


Complex 


cosine 












Square root 


(Arg)^/2 




CSQRT 


Complex 


Complex 



Figure 6-2. Mathematical Subroutines 
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THE CALL STATEMENT 



Mathematical Subroutines 



The CALL statement is used to refer to a SUBROU- 
TINE subprogram. 



General Form 



CALL SUBR ( aj , ag, . . , a^) 

where: 

1. SUBR is the name of a SUBROUTINE subprogram, and 

2. aj, a2, . . , ajj are the n arguments. 



Examples : 

CALLMATMPY (X, 5, 10, Y,7,2) 

CALL QDRTIC (9. 732, Q/4. 536, R-S**2. 0,X1, 

X2) 

The CALL statement transfers control to the sub- 
program and presents it with the actual arguments. 
The arguments may be any of the following: 

1. Any type of constant. 

2. Any type of subscripted or non -subscripted 
variable . 

3. An arithmetic or a logical expression. 

4. Alphameric characters. Such arguments 
must be preceded by nil, where n is the count 
of characters included in the argument, e.g. , 
9HEND POINT. Note that blank spaces and 
special characters are considered in the 
character count when used in alphameric 
fields . 

5. The name of a FUNCTION or SUBROUTINE 
subprogram. 

The arguments presented by the CALL statement 
must agree in number, order, type, and array size 
(except as explained under the DIMENSION state- 
ment) with the corresponding arguments in the SUB- 
ROUTINE statement of the called subprogram. 



The mathematical subroutines are FUNCTION sub- 
programs and are listed in Figure 6-2. 

Machine Indicator Tests 

In the following list of machine indicator test sub- 
routines, assume that i is an integer expression 
and that j is an integer variable. These subroulines 
are referenced by CALL statements. 

SLITE (i): If i = 0, all sense lights will be turned 
off. If i =1,2,3, or 4, the corresponding sense 
light will be turned on. 

SLITET (i, j): Sense light i will be tested and turn- 
ed off. The variable j will be set to 1 if i was on, 
or j will be set to 2 if i was off. 

SSWTCH (i, j): Sense switch i is tested and j is set 
to 1 if i was down and j is set to 2 if i was up. 

OVERFL (j): j is set to 1 if a floating point overflow 
condition exists, or j is set to 2 if no overflow con- 
dition exists. The machine is left in a no overflow 
condition. 

DVCHK (j) ; If the divide check indicator is on, j is 
set to 1 and the divide check indicator is turned off; 
if the divide check indicator is off, j is set to 2. 

EXIT, DUMP, and PDUMP 

EXIT 

A CALL to the EXIT subprogram terminates the 
execution of any program by returning control to the 
Monitor. 

DUMP 



SUBPROGRAMS PROVIDED BY FORTRAN 

FORTRAN includes several commonly used subrou- 
tines that are available to the programmer. The 
mathematical subroutines that are provided are de- 
fined as FUNCTION subprograms; the subroutines 
provided to test the status of the machine indicators 
(the sense switches and the sense lights) are de- 
fined as SUBROUTINE subprograms. In addition, 
FORTRAN includes the SUBROUTINE subprograms 
EXIT, DUMP, and PDUMP. EXIT terminates job 
execution; DUMP dumps core storage and then ter- 
minates job execution; PDUMP dumps core storage 
and then continues execution. 



A CALL to the DUMP subprogram by the statement 

CALL DUMP (Ai, Bi, F], . . . , An, Bn, Fn) 
causes the indicated limits of core storage to be 
dumped and execution to be terminated by returning 
control to the Monitor. 

1. A and B are variable data names that indi- 
cate the limits of core storage to be dumped; 
either A or B may represent upper or lower 
limits. 

2. Fj^ is an integer indicating the dump format 
desired: 

F = dump in octal 

1 dump as real 

2 dump as integer 

3 dump in octal with mnemonics 
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3. 



4. 



If no arguments are given, all of core storage 
is dumped in octal. 

If the last argument F^ is omitted, it is as- 
sumed to be equal to and the dump will be 
octal. 



PDUMP 

A CALL to the PDUMP subprogram by the statement 

CALL PDUMP (Ai, Bi, Fi, . . . , A^, B„^ F^) 
causes the indicated limits of core storage' to be 
dumped and execution to be continued. The PDUMP 
arguments are the same as the DUMP arguments. 

THE BLOCK DATA SUBPROGRAM 

In order to enter data into a COMMON block during 
compilation, a special subprogam must be written. 
This special subprogram contains only the DATA 
COMMON, DIMENSION, and Type statements assoc- 
iated with the data being defined. 

1. The BLOCK DATA subprogram may not con- 
tain any executable statements. 

2. The first statement of this subprogram must 
be the BLOCK DATA statement. 



General Form 



BLOCK DATA 



All elements of a COMMON block must be 
listed in the COMMON statement even though 
they do not all appear in the DATA statement; 
for example, the variable A in the COMMON 
statement in Figure 6-3 does not appear in 
the DATA statement. 
Data may be entered into more than one 
COMMON block in a single BLOCK DATA 
subprogram. 



BLOCK DATA 

COMMON /ELN/A,B,C/RMG/Y,Z 

DIMENSION B( 4), Z( 3) 

DOUBLE PRECISION Z 

COMPLEX C 

DATA (B(I), 1=1, 4)/l. 1,1.2, 2*1. 3/, C/(2. 4, 3. 769)/, 

Z(1)/7.6498085D0/ 
END 



Figure 6-3 
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CHAPTER 7: THE SPECIFICATION STATEMENTS 



The specification statements perform the following 
functions: they provide information about storage 
allocation and about the constants and variables used 
in the program . 

The DIMENSION Statement 



SUBROUTINE MAYMY (..., R, L, M, ... ) 



General Form 



DIMENSION vi(ii), V2(i2), V3(i3), • • • 
where: 

1. each vm is a subscripted variable, and 

2. each in is composed of 1, 2, or 3 unsigned integer 
constants and/or integer variables (i may be a variable only 
when the DIMENSION statement appears in a subprogram). 



Examples: 

DIMENSION A (10), B(5,15), C(L,M) 

DIMENSION S(10), K(5,5,5), G(IOO) 
The DIMENSION statement provides the informa- 
tion necessary to allocate storage for arrays in the 
object program. The DIMENSION statement defines 
the maximum size of arrays. 

1. Each variable that appears in subscripted 
form in the source program must appear in a 
DIMENSION statement contained in the source 
program; however, if the dimension informa- 
tion for a variable is included in a COMMON 
statement in the source program, it must not 
be included in a DIMENSION statement. 

2. The DIMENSION statement must precede the 
first appearance of each subscripted variable 
in an executable or DATA statement for which 
it specified the size. 

3. A single DIMENSION statement may specify 
the dimensions of any number of arrays. 

4. Dimensions specified in a COMMON state- 
ment are subject to all the rules for the 
DIMENSION statement. 

Adjustable Dimensions 

The name of an array and the constants that are its 
dimensions may be passed as arguments in a sub- 
program call. In this way, a subprogram may per- 
form calculations on arrays whose size is not deter- 
mined until the subprogram is called. Figure 7-1 
illustrates the use of adjustable dimensions. 

1 Variables may be used as dimensions of an 
array only in the DIMENSION statement of a 
FUNCTION or SUBROUTINE subprogram. 
For any such array, the array name and all 
the variables used as dimensions must appear 
as arguments in the FUNCTION or SUBROU- 
TINE statement. 



DIMENSION..., R(L,M),, 



DO 100 1=1, L 



Figure 7-1 

2. The adjustable dimensions may not be altered 
within the subprogram. 

3. The absolute dimensions must be specified 

in a DIMENSION statement of the calling pro- 
gram. 

4. The calling program passes the specific di- 
mensions to the subprogram. These specific 
dimensions are those that appear in the DI- 
MENSION statement of the calling program. 
Variable dimension size may be passed 
through more than one level of subprogram. 

The COMMON Statement 



General Form 



COMMON a, b, c, . . . hi A, e, f , . . . /s/g, h, . . . 

where: 

1 . a, b, ... are variables that may be dimensioned, 

and 

2. hi , /s/, . . . are variables that are block names. 



Examples: 

COMMON A,B,C/X/Q,R/YY/M,P,Q 

COMMON /Z/G, H, J//D, F 

Variables, including array names, appearing in a 
COMMON statement are assigned locations relative 
to the beginning of a particular common block. This 
COMMON area may be shared by a program and its 
subprograms. 

1. If the variables appearing in a COMMON state- 
ment contain dimension information, they must not 
appear in a DIMENSION statement. 

2. The locations in the COMMON area are as- 
signed in the sequence in which the variables appear 
in the COMMON statement, beginning with the first 
COMMON statement of the program. 

3. Elements placed in COMMON may be placed 
in separate blocks. These separate blocks may 
share space in core storage at object time. Blocks 
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are given names and those with the same name oc- 
cupy the same space. 

4. COMMON Block Names. The symbolic name 
of a block, which is 1-6 alphameric characters the 
first of which is alphabetic, precedes the variable 
names belonging to the block. The block name is 
always embedded in slashes, e.g. , /BB/. It must 
not be the same as the name of any other subroutine 
which is part of the same job. There are two types 
of COMMON blocks: blank and labeled. 

a. Blank COMMON is indicated either by omit- 
ting the block name if it appears at the be- 
ginning of the COMMON statement, or by 
preceding the blank COMMON variable by 
two consecutive slashes. 

b. Labeled COMMON is indicated by preceding 
the labeled COMMON variables by the block 
name embedded in slashes. 

5. The field of entries pertaining to a block name 
ends with a new block name or with the end of the 
COMMON statement or with a blank COMMON des- 
ignation. 

6. Block name entries are cumulative throughout 
the program. For example, the COMMON state- 
ments 

COMMON A,B,C /R/D,E/S/F 
COMMON G , H/R/I/S/P 
have the same effect as the statement 

COMMON A, B, C, G, H/R/D, E, I/S/F, P 

7. Blank COMMON may be any length. Labeled 
COMMON must conform to the following size re- 
quirement: all COMMON blocks of a given name 
must have the same length in all the programs that 
are executed together. 

8. Variables brought into a COMMON block 
through EQUIVALENCE statements may increase 
the size of the block (see the section, "The Equiva- 
lence Statement"). 

9. Two variables in COMMON may not be made 
equivalent to each other, directly or indirectly. 

10, A double-word variable in COMMON must 
be placed such that its high-order part is an even 
number of words away from the first element in 
COMMON. 

The EQUIVALENCE Statement 



General Form 



EQUIVALENCE (a, b, c, ...), (d, e, f, ...),.. . 
where: 

a, b, c, d, e, f , . . . are variables that may be singly sub- 
scripted; this subscript must be an integer constant. 



Example: 

EQUIVALENCE (A,B(1), C(5)), (D(17), E(3)) 



The EQUIVALENCE statement controls the alloca- 
tion of data storage by causing two or more varia- 
bles to share the same core storage location. 

1. An EQUIVALENCE statement may be placed 
anywhere in the source program. Each pair of 
parentheses in the statement list encloses the names 
of two or more variables that are to be stored in 
the same locations during execution of the object 
program; any number of equivalences (i.e. , sets of 
parentheses) may be given. 

2. In an EQUIVALENCE statement, the meaning 
of C(5) is "the fourth storage location following the 
one that contains C. " In general, C(p) is defined 
for p > to mean the (p-l)th location after the begin- 
ning of the C array; i.e. , the pth location in the 
array. If p is not specified, it is taken to be 1. 

Thus, in the above example, the EQUIVALENCE 
statement indicates that A , and the B and C arrays 
are to be assigned storage locations so that the ele- 
ments A, B(l), and C(5) are to occupy the same loca- 
tion. In addition, it specifies that D(17) and E(3) are 
to share the same location. 

3. Quantities or arrays that are not mentioned in 
an EQUIVALENCE statement will be assigned unique 
locations. 

4. Locations can be shared only among variables, 
not among constants. 

5. The sharing of storage locations requires a 
knowledge of which FORTRAN statements will cause 
a new value to be stored in a location. There are 
four such statements: 

a. Execution of an arithmetic statement 
stores a new value in the variable on the 
left side of the equal sign. 

b. Execution of an ASSIGN i TO n statement 
stores a new value in n. 

c. Execution of a DO statement or an implied 
DO in an I/O list will store a new indexing 
value. 

d. Execution of a READ statement stores 
new values in the variables mentioned in 
the input list. 

6. Variables brought into a COMMON block 
through EQUIVALENCE statements may increase 
the size of the block indicated by the COMMON 
statements, as in the following example: 

COMMON /X/ A, B, C 

DIMENSION D(3) 

EQUIVALENCE (B, D) 
The layout of core storage indicated by this ex- 
ample (extending from the lowest location of the 
block to the highest location of the block) is: 

A 

B, D(l) 

C, D(2) 
D(3) 
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7. Since arrays must be stored in consecutive 
forward locations, a variable may not be made equi- 
valent to an element of an array in such a way as to 
cause the array to extend beyond the beginning of 
the COMMON block. For example, the following 
coding is invalid: 

COMMON /X/A, B, C 
DIMENSION D(3) 
EQUIVALENCE (B, D(3)) 
because it would force D(l) to precede A, as fol- 
lows: 

D(l) 
D(2) 
D(3) 



A, 
B, 
C 



8. The rule for equivalencing double-word vari- 
ables with single-word variables is: 

In COMMON, the effect of the EQUIVALENCE 
statements must be such that the high-order word of 
any double->-word variable is an even number of loca- 
tions away from the start of the COMMON block. 

In non-COMMON, the effect of the EQUIVALENCE 
statements must be such that the high-order word 
of any double-word variable is an even number of 
words away from the start of any other double- 
word variable linked to it through EQUIVALENCE 
statements. 

9. Two variables in one COMMON block or in 
two different COMMON blocks may not be made 
equivalent. 

The Type Statements 

The type of a variable or function may be specified 
by means of one of the six Type statements. 



General Form 



INTEGER a,b,c,... 

REAL a, b,c, . . . 

DOUBLE PRECISION a, b, c, . . . 

COMPLEX a, b, c, . . . 

LOGICAL a,b, c, ... 

EXTERNAL a, b, c, . . . 

where: 

a, b, c, . . . are variable or function names appearing within 
the program. 



Examples: 

INTEGER BIXF, X, QF, LSL 
REALIMIN, LOG, GRN, KLW 
DOUBLE PRECISION Q, J, DS IN 
COMPLEX DZ,DX,M,T 
LOGICAL F, G, LWG 
EXTERNAL SIN, MATMPY, INVTRY 



The variable or function names following the type 
(INTEGER, REAL, etc.) in the Type statement are 
defined to be of that type, and remain that type 
throughout the program; the type may not be changed. 

Note that LSL and GRN need not appear in their 
respective type statements, since their type is 
implied by their first characters. 

1. The appearance of a name in any Type state- 
ment, except EXTERNAL, overrides the implicit 
type assignment. 

2. Variables that are subprograms names that 
are arguments of other subprograms must appear 
in EXTERNAL statements. 

3. A name may appear in two Type statements 
only if one of them is EXTERNAL. 

4. The Type statement must precede the first 
use of a name in any executal)le statement or DATA 
statement in the program. 

5. A name declared to be of a given type may 
assume only the values of a constant of the same type. 

The DATA Statement 

Data may be compiled into the object program by 
means of the DATA statement. 



General Form 



DATA list/di , d2, . . . , dn/ , list/di , 62, k*d3, . . . , dm/ , . 
where: 

1. list contains the names of the Ntariables being 
defined, 

2. d is the information literal, and 

3. k is an integer constant. 



Examples: 

DATA R, Q/14. 2, 3HEND/, Z/O777777700001/ 
DATA (B(I),C(I), 1=1, 40, 2)/2. 0,3. 0,38*100.0/ 

1. List. Subscripted variables may appear in 
the Hst of values. Where a subscript symbol is 
used, it must be under control of DO- implying 
parentheses and associated parameters. Subscripts 
not so controlled must be integer constants. The 
DO-defining parameters must be integer constants. 

2. k. k may appear before a d-field to indicate 
that the field is to be repeated k times, k must be 
followed by an * to separate it from the field to be 
repeated. 

3. d. The d-literals may take any of the three 
following forms: 

a. Integer, real, double -precision, and com- 
plex constants. 

b. Alphameric characters. The alphameric 
field is written as nH followed by n alpha- 
meric characters. Each group of six alpha- 
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c. 



4. 



meric characters form a word. If n is not 
a multiple of six, the remaining characters 
are left -adjusted in the word and the word 
is filled out with BCD blanks. Blanks are 
significant in alphameric fields. 
Octal digits. The octal field is written O 
followed by 1 - 12 signed or unsigned octal 
digits. 

There must be a one-to-one correspondence 



between the list items and the data literals. For 



example, if it is desired to define 16 alphameric 
characters, say 16HDATAbTObBEbREADb, as a 
variable, say G, then G must be dimensioned to 
cover at least 3 locations, 

5. The BLOCK DATA subprogram (see page 27), 
which includes a DATA statement, compiles data 
into the common area of the program. 

6. When DATA defined variables are redefined 
during execution, these variables will assume their 
new values regardless of the DATA statement. 
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APPENDIXES 



APPENDIX A: SOURCE PROGRAM STATEMENTS 
AND SEQUENCING 

The order in which the source program statements 
of a FORTRAN program are executed follows these 
rules: 

1. Control originates at the first executable 
statement. The specification statements 
and the FORMAT, FUNCTION, SUBROU- 
TINE, and BLOCK DATA statement, are 
non-executable; in questions of sequencing, 
they can be ignored. 

2. If control was with statement S, then control 
will pass to the statement indicated by the 
normal sequencing properties of S. If, how- 
ever, S is the last statement in the range of 
one or more DOs that are not yet satisfied, 
then the normal sequencing of S is ignored 
and DO-sequencing occurs. 

The normal sequencing properties of each FOR- 
TRAN statement follows . 

Statement Normal Sequencing 



a = b 

ASSIGN i TO n 
BACKSPACE i 
BLOCK DATA 
CALL 

COMMON 
COMPLEX 
CONTINUE 



Next executable statement 

Next executable statement 

Next executable statement 

Non-executable 

First statement of called 

subprogram 

Non- executable 

Non-executable 

Next executable statement 



Statement 


Normal Sequencing 


DATA 


Non-executable 


DIMENSION 


Non-executable 


DO 


DO-sequencing, then the next 




executable statement 


DOUBLE PRECISION 


Non-executable 


END 


Terminates program 


END FILE 


Next executable statement 


EQUIVALENCE 


Non-executable 


EXTERNAL 


Non-executable 


FORMAT 


Non-executable 


FUNCTION 


Non-executable 


GO TO n 


Statement n 


GO TO n, (ni, n2, . . . , nm) 


Statement last assigned to n 


GOTO(ni,n2,...,njn),i 


Statement nj 


IF (a) ni,n2,n3 


Statement ^ifT^2*^3 *^ * "^ *^> 




a=0, or if a > 0, respectively 


IF (t) s 


Statement s or next executable 




statement if t is true or false, 




respectively 


INTEGER 


Non-executable 


LOGICAL 


Non-executable 


PAUSE 


Next executable statement 


PRINT 


Next executable statement 


PUNCH 


Next executable statement 


READ 


Next executable statement 


REAL 


Non-executable 


RETURN 


The first statement, or part of a 




statement, following the reference 




to this program 


REWIND 


Next executable statement 


STOP 


Terminates the program 


SUBROUTINE 


Non- executable 


WRITE 


Next executable statement 
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